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SUMMARY 

This  report  documents  the  development  of  a  mathematical  model  ^hich 
describes  the  gunner's  tracking  performance  in  an  AAA  tracer-directed  manual 
firing  task.  Reduced-order  observer  theory  is  applied  to  design  this  model 
for  the  underlying  linear  time-varying  antiaircraft  artillery  system.  The 
model  consists  of  a  reduced-order  observer,  a  linear  feedback  controller, 
and  a  stochastic  remnant  element.  Both  the  tracking  and  the  tracer  errors 
are  considered  measurable.  The  tracer  dynamics  enter  the  system  in  a 
delayed  fashion.  The  Average  Approximation  Method  is  used  to  solve  the 
closed-loop  delay  differential  equations.  A  least-squares  minimization 
algorithm  is  derived  to  identify  the  model  parameters  systematically. 
Simulation  results  on  model  predictions  versus  empirical  data  over  several 
input  trajectories  are  included.  The  gunner  model  can  adequately  describe 
human  response  in  this  compensatory  tracking  and  firing  task.  The  gunner 
model  developed  here  can  be  used  in  larger  attrition  models  to  evaluate  the 
survivability  of  aircraft  in  tactical  engagement  scenarios. 
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Section  I 
INTRODUCTION 


The  tracking  performance  of  a  human  operator  in  an  antiaircraft  artillery 
(AAA)  system  has  been  studied  for  several  years.  Several  approaches  have 
been  proposed  and  proven  successful  in  modeling  human  response  in  compen¬ 
satory  tracking  tasks.  Among  these  are  the  McRuer  Crossover  Model  (McRuer 
and  Krendel,  1974),  Optimal  Control  Model  (Kleinman  et  al.,  1971,  1974),  PID 
Structure  Modified  Optimal  Control  Model  (Phatak  et  al.,  1976),  and  Observer 
Model  (Kou  et  al.,  1978).  The  specific  gunnery  task  modeled  involves  direct 
manual  rate  control  of  the  gun  turret  with  the  gunner  using  an  optical 
sighting  system  having  a  line  of  sight  coincident  with  the  gun  pointing 
angle.  In  this  system  mode,  radar  is  not  used  for  either  azimuth,  eleva¬ 
tion,  or  range  data.  Rather,  the  gunner,  using  visual  feedback  from  the 
antiaircraft  artillery  tracer  rounds  he  has  fired,  continuously  adjusts 
weapon  pointing  in  azimuth  and  elevation  to  minimize  the  tracer-to-target 
error;  hence,  increases  his  chances  of  putting  "hits"  on  the  target.  There¬ 
fore,  the  subject  plays  both  the  role  of  a  tracker  and  a  lead  angle  computer. 
This  tracking  task  is  greatly  complicated  by  the  inclusion  of  lead  angle 
estimation.  An  additional  measurement  channel  is  available  to  the  gunner  to 
supplement  the  task  through  the  miss  distance  of  tracer  rounds  from  the 
target.  The  tracer  dynamics  enter  the  man-machine  system  in  a  delayed 
fashion  and  complicate  the  development  of  a  simple  but  faithful  human  opera¬ 
tor  model.  In  this  report,  Luenberger  observer  theory  is  applied  to  design 
an  antiaircraft  gunner  model  that  consists  of  a  reduced-order  observer  and  a 
controller.  A  remnant  function  that  lumps  all  the  random  effects  from 
measurement  noise  and  human  neuroraotor  response  noise  is  assumed  to  be 
Gaussian  with  its  covariance  being  a  function  of  target  velocity  and  accel¬ 
eration.  The  manned  AAA  simulation  experiment  was  conducted  at  the  Air 
Force  Aerospace  Medical  Research  Laboratory,  Wright-Patterson  AFB,  Ohio. 

The  empirical  tracking  and  tracer  data  were  used  to  identify  parameters  of 


empirical  data  over  various  input  target  trajectories. 
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Section  II 

AAA  TRACKING  AND  FIRING  SYSTEM 

In  a  manual  tracking  and  firing  task,  the  human  operator  (gunner)  perceives 
both  the  tracking  error  and  the  tracer  error  on  a  two-dimensional  visual 
display.  The  tracking  error*  e^  is  the  difference  between  the  target  angle 
0T  and  the  barrel  angle  0g,  while  the  tracer  error  e2  is  the  difference 
between  the  target  angle  6p  and  the  projectile  ending  angle  6p.  For  the 
simulation  the  projectile  flight  path  ended  at  the  range  of  the  target; 
each  individual  tracer  round  was  blanked  at  this  point  6p  and  disappeared 
from  the  display.  This  is  depicted  in  Figure  1.  The  two  axes  on  the  dis¬ 
play  represent  the  azimuth  and  elevation  components  of  the  tracking  and  the 
tracer  errors,  respectively.  The  gunner's  objective  is  to  align  the  pro¬ 
jectile  ending  angle  0p  to  the  target  angle  0  .  or  to  minimize  the  magni¬ 
tude  of  the  tracer  error  to  simulate  a  close  hit  on  the  target. 


The  gunner  performs  both  the  tracking  and  firing  task  along  both  the  azimuth 
and  elevation  axes.  Since  the  structure  of  the  underlying  systems  for 
azimuth  and  elevation  axes  is  similar,  they  can  be  decomposed  and  treated 
separately.  The  only  difference  comes  from  the  azimuth  tracking  system 
being  nonlinear  due  to  a  visual  correction  factor  cos(0D)„T  .  Figure  2  is 
the  block  diagram  of  an  AAA  tracking  and  firing  system. 

ELEVATION 

The  elevation  barrel  dynamics  and  rate  control  system  can  be  represented  by 
the  transfer  function  T-^g 


,  v  e!B(s)  1.34  (s+16.875) 

ilB^S^  U  (s')  2 

lvs;  s(0. 926s  +S+16. 875) 


(1) 


Based  o-i  a  frequency  domain  analysis  of  the  input  target  trajectories  (Kou 
and  Glass,  1977),  it  was  found  that  the  frequency  bandwidth  of  all  the 


*Also  referred  to  as  the  "lag"  angle  later  in  this  report. 


Figure  1.  Visual  Display  of  Tracking  and  Tracer  Errors 


trajectories  used  in  tnis  study  is  approximately  0.2  Hz.  Hence,  the 
above-mentioned  transfer  function  can  be  approximated  by 


1.34 

s 


(2) 


to  simplify  the  model  design  and  subsequent  analysis.  The  corresponding 
input-output  relation  of  the  elevation  barrel  system  becomes 


^1B(t)  =  l-34  ^1(t)  (3) 

The  elevation  projectile  ballistics  can  be  described  by  the  following 

equation  (see  Appendix  A):  ,j 

| 

2  $ 
O'ip(t)  =  01B(t-T)  -  0.001  (5. 2t+0.486t  )  cos  0X  (t-r)  (4)  J 

t 

■4 

y 

i 

where  t  is  the  time-of-flight  of  the  projectile.  The  second  term  approxi¬ 
mates  the  elevation  drop  of  a  projectile  due  to  gravitation.  The  t  is  deter¬ 
mined  by  the  target  range  R(t)  (see  Appendix  A) 


i 

5 

\ 


R(t) 

930  -  0.19  R(t) 


for  R(t)  <  2877  meters 


and  thus  is  a  time-varying  quantity. 


(5) 


AZIMUTH 


The  azimuth  barrel  dynamics  and  rate  control  system  can  be  represented  by 
the  following  transfer  function 


v2B(s)  =  1.28  (s+5.312) 

U  (s')  ? 

2K  J  s(0.459s“+s+5.312' 


(6) 


Similar  to  the  elevation  case,  the  azimuth  transfer  function  can  be  approxi¬ 
mated  by 


t2B(s) 


1.28 


(7) 
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The  corresponding  input-output  relation  of  the  azimuth  barrel  plant  becomes 


°2B  (t)  =  1,28  U2(t) 


Since  there  is  no  gravitational  drop  in  the  azimuth  axis,  the  azimuth  pro¬ 
jectile  ballistics  can  be  described  simply  by 


e2P  (t)  =  °2B  (t“T) 


The  state  vector  is  introduced  here 


=  [yo.  x.2(t),  x.3(t)jT, 


"T"  means  "the  transpose  of,"  with 


x.,(t)  =  q4M  -  e._(t). 


x .At)  =  e  (t)  -  eip(t) 


xi3(t)  =  o.T(t), 


where 


i  =  1,  2* 


*If  not  otherwise  specified,  the  first  subscript  index  i  represents  the 
elevation  (i  =  1)  or  azimuth  axis  (i  =  2)  while  the  second  index  rep 
sents  the  i-th  element  or  row  of  a  matrix. 


11 


The  following  system  and  measurement  equations  can  be  derived  from 
Equations  (3)  —  (4) ,  (8)  —  (9 ) . 


with 

bx  =  -1.34 
b2  =  -1.28 

C1  =  1 

c2  =  C0S  9lB(t) 

ex(t)  =  -1.34  x  (1— x)  x  jl  +  (0.0052t+0.000486t2)  sin  0lB(t-T)] 
e2(t)  =  -1.28  x  (1-t) 

gl(t)  «  (0. 0052+0. 000972T)  x  T  x  cos  91B(t-x) 
g2(t)  =  0 
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8f-»  : 


0  ,  u  ,  y^,  and  y ^  denote  the  elevation  or  azimuth  components  of  the 

target  acceleration,  the  gunner’s  control  output  and  the  observed  tracking 
error  (lag  angle)  and  tracer  error,  respectively.  By  definition,  0  (t) 

can  be  expressed  by  ©.^(t)  -  x^(t).  Thus,  the  underlying  AAA  tracking  and 
firing  system  is  not  only  a  time-varying  system,  but  also  a  nonlinear 
system  with  a  time-varying  delay.  In  order  to  have  the  problem  well-posed, 
several  approximations  are  made  to  simplify  the  system  equations. 

Note  that  0,  (t)  is  not  directly  measurable  and  T  is  less  than  or  equal  to 
IB 

7.5  seconds  in  the  simulation  experiment,  and  we  may  neglect  g^(t)  and 
approximate  e^(t)  by 

e  (t)  =  -1.34  (l-x)ei  (12) 

where  8^  is  a  parameter  to  be  identified  with  other  model  parameters. 

If  we  introduce  a  transformation  on  the  states  x^  and  x^  by  xlj  = 

x'  =  c.  x.„  then  Equations  (10)-(11)  can  be  rewritten  as  follows.  Let  8* 

i2  l  i2  *■ 

=  1, 


x!  =  A!  x!  +  B!  u.(t)  +E!  u.  (t-T)  +  F  0.T(t)  (13) 

1  11  XI  11  1  11 

y^(t)  =  £1^  x'  i(t)  (14) 


wiiere 


a!  = 

— i 

-i 

c  .c . 

1  1 

0 

0 

‘  -1 
cici 

ci 

c . 

l 

5i- 

cibi 

0 

0 

0 

0 

0 

c.b.  (l-i)g. 
11  1 


4  =  0 


i  =  1,2 


Notice  that  g^(t)  is  dropped  from  Equation  (10)  and  e^(t)  is  approximated 
by  the  expression  shown  in  Equation  (12).  The  simplified  system  equations 
represent  a  linear  time-varying  system  with  a  time-varying  delay  in  the 
control. 


Section  III 

AN  AAA  OBSERVER  GUNNER  MODEL 


The  function  of  the  gunner  in  an  AAA  system  can  be  regarded  as  consisting 
of  two  stages  to  be  modeled.  In  the  first  stage,  the  gunner  observes  con¬ 
tinuous  signals  and  makes  an  estimate  of  system  states  (e.g.,  target 
velocity)  based  on  his  internal  model  of  target  motion.  In  the  second 
stage,  the  gunner  uses  the  observed  and  estimated  states  to  form  and  exercise 
a  proper  control  action  so  that  a  certain  objective  can  be  achieved.  The 
former  one  corresponds  to  an  estimation  process,  while  the  latter  corresponds 
to  a  control  process.  The  reduced-order  observer  is  used  in  conjunction 
with  a  linear  state  variable  feedback  (l.s.v.f.)  control  law  to  model  the 
gunner's  function.  The  structure  of  an  observer  gunner  model  is  shown  in 
Figure  3. 


Since  the  gunner  does  not  have  full  knowledge  of  the  target  dynamics,  we  may 
assume  that  the  gunner  perceives  the  target  motion  (or  internal  model)  as  a 
constant  velocity  process,  i.e.,  @T  =  0.  Therefore,  the  equations  repre¬ 
senting  the  gunner's  internal  model  of  the  tracking  and  firing  system  can  be 

described  by 


x!  =  A!x!  +  B!  u.  (t)  +  E ’  u  (t-T) 

—1  — 1—1  “I  1 


(15) 


*i 


C'.x , 


i  =  i,: 


(16) 


Note  that  C!  is  a  constant  matrix  and  both  x^,  xi2  are  measurable,  the  only 


state  that  needs  to  be  estimated  in  order  to  implement  an  l.s.v.f.  control 
law  is  x.  ,  the  target  velocity.  Reduced-order  observer  (Luenberger,  1971) 
is  applied  to  generate  an  estimate  xi3  of  x„3  by  using  Equations  (15)  and 
(16).  The  estimate  x±3  satisfies  the  following  equation: 
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i 


a 


HUMAN  OPERATOR  MODEL  REMNANT 


Xi3  ~  (kil+ki2^  Xi3^  +  kilyil^  +  ki2yi2^ 


12J12 


'  kilCiCi_1  yil(t;  -  ki2CiCi"lyi2(t) 


-  B±  (1-T)  ki2b±  u±(t-T)  -  kilbiui(t)  i  =  1,2 


where  and  k^  are  observer  gains. 

As  mentioned  earlier  in  this  section,  the  second  function  of  the  gunner  is 
to  form  a  control  action  based  on  the  observed  tracking  and  tracer  error  as 
well  as  the  estimated  target  velocity.  The  gunner  uses  a  control  manipulator 
called  an  "H-grip"  to  position  the  gun  barrel/optical  sight  in  azimuth  and 
elevation;  this  rate  control  system  is  continuously  used  by  the  gunner  to 
minimize  the  error  between  the  tracer  rounds  and  the  target  (i.e.,  minimize 
*2  of  Figure  3).  The  objective  of  the  gunner  is  to  minimize  the  tracer 
error  so  that  a  maximum  probability  of  hit  could  result.  In  other  words, 
the  gunner's  response  in  the  second  stage  would  be  to  stabilize  the  under¬ 
lying  system,  especially  the  tracer  error  xi2(t).  A  natural  control  func¬ 
tion  which  may  achieve  this  objective  would  be  a  linear  feedback  control 
law.  Let  us  consider  an  l.s.v.f.  control  of  the  form 


u.(t)  =  F.  x!  (t)  +  v.(t) 
l  —l—i  l 


i  =  1,2 


where  T,  is  a  vector  of  controller  gains  to  be  identified,  x!  is  a  vector 

—l  i 

of  measurable  states  and  estimated  state. 
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A  remnant  function  v^(t)  is  introduced  here  to  take  into  account  the  random 
effects  due  to  the  measurement  noise,  the  neuromotor  noise,  target  uncer¬ 
tainty,  modeling  error,  difficulty  of  tracking,  etc.  It  is  modeled  as  a 
white  noise  with  zero  mean  and  with  covariance  function 


E  v.  (t)  v i  (s)  =  ot^  +  a±2  0iT  (t)  +  ai3  0iT(t)  5(t-s) 


i  =  1,2 


fo-  all  t  and  s.  ol  are  nonnegative  model  parameters  to  be  determined. 


6^^  and  0  are  estimated  target  angle  rate  and  acceleration,  respectively. 


The  model  Equations  (17)  and  (18)  describe  the  gunner's  response  in  the 
estimation  and  the  control  stage,  respectively.  The  overall  equations  of 
the  man-in- the-loop  AAA  tracking  and  firing  system  are  obtained  by  com¬ 
bining  Equations  (17)  and  (18)  with  Equations  (13)  and  (14)  of  the  actual 
tracking  and  firing  system.  If  we  define  a  new  state  vector 


X. 

— l 


x 


i3 


then  the  state  equation  of  the  closed-loop  system  becomes 


(20) 


X. (t)  +  D.  X. (t-T)  +  F.0 
— i  —l—i  —l 


iT 


+  B.v.(t)  +  E.v.(t-T) 
—li  —li 


(21) 


where 


Section  IV 

A  LEAST-SQUARES  IDENTIFICATION  PROGRAM 


Since  the  last  two  components  of  are  independent  from  the  first  two 
components  as  indicated  in  Equation  (21),  the  system  equations  can  be 
decoupled  into  the  following  form: 


Ki3(t)  -  e.T<t) 


x.4(t)  -  -k^  *.4<t>  +  eJT<t) 


x . ( t )  =  (c.c.  '  +  b.c.Y.,)  x.-,(t)  +  b.c.  Y.0x..(t) 
ll  11  1  1  ll  ll  l  l  i2  i2 


+  b.c.v.(t)  +  f.,  (t) 
ill  ll 


xi2(t)  ~ 


’  Vi"1  xi2(t)  +  bici6i<1-T)  [y^jCt-T) 

+  Yi2xi2^t-T^  +  Vi(t-T)]+  fi2(c^ 


£il(t)  *  <1+biYi3>  Ci  Xi3(t)  -  biCiYi3Xi4(C) 

£12(t)  ’  cix13<t)  +  bi=iBi  (1-T)  Yi3[xi3(t-T)  -  Xi4(£-T)]  (24) 


i  -  1,2 


Note  that  Equation  (23)  represents  a  linear-time  varying  system  with  time- 
varying  delay.  Since  an  explicit  solution  for  this  type  of  delay  dif¬ 
ferential  equation  is  usually  not  available,  the  Average  Approximation 
Method  is  applied  to  approximate  Equation  (23)  by  the  following  ordinary 
differential  equation  (Banks  and  Burns,  1978). 

•  m  N  N  n 

wj(t)  =  AT (t)  W*(t)  +  M±  ^(t) 


(25) 


&, 


2  (N+l) 


0 

1 


0 

0 


fu(t)  +  bjCj.^t) 


fi2<t)  +  b.c.e.d-x)  vi(t-r) 


i  =  1,2 


and 

-i  =  Xil(t)’  xi2(t)»  xil(t_  IT  T)’  Xi2  (t-  “IT  T)’  *  *  ‘ 

Xii(C-T),  x.2(t-T)jT 

is  a  2(N+1)  x  1  vector  consisted  of  the  original  states  x^(t),  xi2(t)  and 
the  approximated  delayed  states  of  order  N. 


Based  on  a  comparison  study,  the  time  history  of  x^(t)  and  x^2(t)  fell 
within  5  percent  of  each  other  when  N  =  1  was  used  versus  N  =  10.  Therefore 
N  =  1  was  adopted  to  save  computer  memory  and  reduce  computation  time.*  For 
N  =  1 ,  we  have 


wj(t)  =  A*(t)  wj(t)  +  £.(t) 


(26) 


^Memory  saving  is  up  to  50  percent  and  computation  time  can  be  reduced  up 
to  75  percent. 
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Vi  +  Wil  biVi2 


Vi  btciYil8i(1-T)  bicl'i2eia“,) 


In  order  to  identify  model  parameters,  we  first  take  the  expectation  of 
Equation  (26)  to  obtain  the  equation  governing  the  mean  of  states 


(t)  =  A*(t)  W^t)  +  mJ  f^t) 


where 


W^(t)  =  E  jx11<t)j  ,  E|xi2(t)j,  ^(t-x)],  E|xi2(t-T)j 


^(t)  =  f±1(t),  fi2(t) 


Notice  that  the  remnant  terms  v. (t)  and  v. (t-x)  disappear  because  of  the 

11  —1 

assumption  of  its  zero  mean.  The  first  and  second  component  of  W^,  w^  and 
wi2’  are  t*ie  mot*e^  Prediction  of  ensembled  mean  of  tracking  and  tracer 
error,  respectively. 


On  the  other  hand,  the  covariance  matrix  P^(t)  satisfies  the  following 
equation: 


P.(t)  =  A*(t)  P.(t)  +  P.(t)  A ( t )  +  D.(t)  £.(t)  D.T(t)  (28) 


where 


P.(t)  =  E 


wj(t)  -  W^( t) 


wj(t)  -  w^Ct) 


D.(t}  = 


b  .c . 

1  1 

0 


0 

Lo 


b.c.B.(l-r) 

ill 

0 

c 


a..  +  a. 

ll  x2 


oiT(t)  ;  +  ai3 


0  ail  +  ai2 


OiT(t-T) 


oiT(t) 


+ai3  OiT(t-T) 


The  first  diagonal  element  of  I\(t)  is  the  square  of  the  ensembled 

standard  deviation  of  tracking  errors.  Similarly,  the  second  diagonal 

S'  #  _ 

element  p  „„  of  P  (t)  is  the  square  of  the  ensembled  standard  deviation  oi 

i22  —i 

tracer  errors.  By  solving  the  matrix  differential  Equation  (28),  and 
defining 


1  1 

sil(t)  =  Pill(t)i’  Si2(t)  =  pi22(t)2  (29) 

the  model  prediction  of  standard  deviation  of  tracking  and  tracer  errors  is 
obtained. 

The  parameters  of  the  gunner  model  are  determined  via  a  least-squares  curve¬ 
fitting  identification  program.  The  reference  curves  to  be  fitted  in 
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the  curve-fitting  program  are  obtained  from  empirical  tracking  and  tracer 
data  collected  in  the  manned  simulation  experiments.  These  experiments  were 
conducted  on  an  AAA  simulator  at  the  Air  Force  Aerospace  Medical  Research 
Laboratory.  Several  simulated  nonmaneuvering  (i.e.,  flyby)  and  maneuvering 
(i.e.,  recon)  aircraft  trajectories  were  used  as  target  trajectories  for 
these  experiments.  Let  x^(t)  and  x^(t)  be  the  empirical  means  of  tracking 
and  tracer  errors,  s^(t)  and  s^ (t)  be  the  corresponding  standard  devia¬ 
tions.  These  empirical  means  and  standard  deviations  were  obtained  by 
averaging  and  computing  the  variance  of  the  empirical  data  of  forty  experi¬ 
mental  simulation  runs  with  the  same  target  trajectory  and  the  same  subject. 


The  parameters  were  identified  by  minimizing  the  cost  function  J(k,  JT, 
a_,  6,)  defined  as  follows: 


min  J^(k,  6,  f_,  a) 


k,8,£,a 


« 


^(t) 


(30) 


i  =  1,2 

where  t  is  the  time  when  the  first  tracer  round  reaches  the  range  of  the 
o 

target,  tf  is  the  time  when  the  last  tracer  round  is  fired.  The  parameters 
to  be  determined  are  V,  8,  _T  and  a  as  defined  in  Equations  (13),  (17)-(19). 
I,  is  a  positive  weighting  factor  chosen  to  be  one  in  the  identification 
runs. 


The  cost  function  consists  of  four  contributing  terms.  These  terms  are, 
respectively: 


•  The  integral  of  the  square  errors  between  the  empirical  mean 
tracking  error  x..^  and  the  model  prediction  ^  ^w^. 


The  square  errors  between  the  empirical  mean  tracer  error  x  0  and 
the  model  prediction  c^  w^- 


The  square  errors  between  the  empirical  standard  deviation  s 

tracking  error  and  the  model  prediction  c.  ^s... 

1  ll 


il 


of 


Finally,  the  square  errors  between  the  empirical  standard  deviation 
si2  of  tracer  error  and  the  model  prediction  c^^s^* 


Therefore,  the  model  parameters  will  be  identified  such  that  model  prediction 
of  tracking  and  tracer  error  and  their  standard  deviations  will  fit  with  the 
empirical  counterparts  simultaneously.  This  is  a  fairly  complex  minimization 
problem  which  would  b?  too  cumbersome  and  inefficient  to  use  the  Gauss- 
Newton  gradient  method.  A  direct  search  method  based  on  Rosenbrock’s  hill¬ 
climbing  minimization  algorithm  (Rosenbrock,  1966)  was  derived  (see 
Appendix  B) .  This  method  adjusts  parameter  values  along  a  set  of  ortho¬ 
normal  axes  iteratively  to  search  for  the  least  cost  J^,  before  a  rotation 
of  axes  is  performed  and  a  new  series  of  searches  are  launched.  The  itera¬ 
tive  process  will  continue  until  the  changes  in  parameter  values  are  smaller 
than  a  prespecified  value.  A  computer  curve-fitting  program  was  written 
which  implements  the  previously  mentioned  minimization  algorithm.  Some 
merits  of  the  identification  program  are  that  it  has  simple  computation 
requirements,  provides  a  systematic  search,  and  gives  reasonably  fast  con¬ 
vergence  speed. 

V 

A  set  of  converged  parameter  values  were  obtained  by  this  identification 
program  for  both  the  elevation  and  azimuth  gunner  model.  These  values  are 
listed  in  Table  1. 


The  preceding  parameter  values  were  obtained  by  applying  the  empirical  data 
from  the  Re^on  flight  path.  Notice  that  the  gunner  model  depends  on  the 
dynamics  of  the  barrel  plant.  For  different  AAA  tracking  systems,  the 
identification  program  can  be  used  to  determine  parameter  values  for  the 
corresponding  gunner  model. 
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Section  V 

COMPUTER  SIMULATION  RESULTS 


The  gunner  model  described  in  the  last  section  was  implemented  on  a  CDC 
CYBER  175  computer  to  simulate  the  man-in-the-loop  AAA  tracking  and  firing 
ta_.k.  The  parameters  of  the  gunner  model  listed  in  Table  1  were  used  in  a 
sifjlation  of  gunner's  response  for  five  typical  operational  trajectories, 
including  both  flyby  as  well  as  maneuvering  trajectories.  In  order  to  avoid 
using  convolution  integration  (to  reduce  computation  time)  in  solving 
Equations  (27)  and  (28),  we  discretize  the  system  equations  into  (for  sim¬ 
plicity  of  notation  we  may  drop  the  subscript  i  and  the  superscript  1  from 
here  on) , 

W=$W+Hf+Lv  (31) 

— n+1  — n  — n  - n  — n  — n  — n 


where 

Hn+i  =  W(cn+l)  with  Cn+1  =  Co  +  (n+1)A 
and  A  =  0.06  seconds. 


and  f  =  f(t  )  as  defined  in  Equation  (24).  Notice  that  v  defined  here  is 
— n  —  n  ~l 

a  random  sequence  with  the  following  properties. 


E[sJ  -  o.  *  [(!„)  <v/]  -  s  a<t„) 

Taking  the  expectation  of  both  sides  of  Equation  (31)  we  obtain 


(32) 


W  _  =  $  W 
— n+1  — n  — n 


H  f 
— n  — n 


(33) 


Denote  P  , ,  as  the  covariance  matrix  of  W  . , ,  then  the  solution  of 
—n+1  -n+1’ 

Equation  (28)  can  be  expressed  in  the  following  discrete  form: 


4+1 


=  $  P  4> 

— n  — n  — n 


+  j  L 
A  — n 


44 


(34) 


The  first  and  second  elements  of  Wn+^  in  Equation  (33)  are  the  model  pre¬ 
diction  of  the  mean  tracking  and  tracer  error  correspondingly.  The  first 
and  second  diagonal  elements  of  P^^  in  Equation  (34)  are  the  model  pre¬ 
diction  of  the  variance  of  tracking  and  tracer  errors. 


I 

< 

t 

k 

;1 
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A  computer  program  was  developed  which  uses  the  recursive  Equations  (33)  and 
(34)  to  simulate  a  closed-loop  AAA  tracking  and  firing  task.  Inputs  to 
the  simulation  program  are  the  time  history  of  range  and  acceleration  of  the 
target  aircraft,  plus  the  initial  angular  ,>sition  and  velocity  of  the 
target.  Outputs  of  the  simulation  program  are  model  predicted  mean  tracking 
error  and  its  standard  deviation.  Five  trajectories  including  flyby  and 
maneuvering  trajectories  were  used  in  this  study  and  are  shown  in  Figure  4a- 
4e.  The  AAA  weapon  system  was  located  at  the  origin  of  the  x-y  plane,  with 
aircraft  altitude  measured  along  the  -z  axis.  The  increment  of  each  of  the 
three  axes  is  1000  feet ,  A  detailed  description  of  the  characteristics  of 
these  trajectories  can  be  found  in  (Rolek,  1977). 


Simulation  results  are  shown  in  Figures  5  through  24  for  the  trajectories 
listed  in  Figure  4.  The  solid  curve  in  these  figures  is  the  empirical  data  ] 

which  are  obtained  by  averaging  the  results  of  40  experimental  simulation  . 

runs.  The  dashed  curve*  is  the  model  prediction  of  ensembled  mean  and 

*Initial  conditions  used  in  these  simulation  runs  are  the  same  as  the  j 

empirical  data.  In  general,  initial  conditions  are  computed  by  empirical  | 

formulas  referred  in  Appendix  A.  | 


'  'Vj*  «1KV**^v^$VV>V 


standard  deviation.  Figures  5  through  8  show  the  results  of  elevation  mean 
and  standard  deviation,  azimuth  mean  and  standard  deviation  of  both  tracking 
errors  and  tracer  errors  for  the  Recon  trajectory.  These  figures  show  that 
the  designed  gunner  model  can  provide  consistent  prediction  of  the  human's 
empirical  tracking  data  as  well  as  the  tracer  error  data  which  result  when 
tracking  error  propagates  through  the  barrel  plant.  The  remaining  figures 
show  that  similar  agreement  between  predicted  and  empirical  data  holds  for 
the  other  four  trajectories,  too.  The  author  concluded  that  for  a  given  AAA 
weapon  system,  the  same  set  of  parameter  values  can  be  used  in  the  gunner 
model  to  predict  the  human  operator's  tracking  and  tracer  errors  for  all 
simulated  target  trajectories.  In  other  words,  the  gunner  model  designed 
here  is  a  predictive  model. 
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igure  13b.  Standard  Deviation  of  Tracking  Error — Elevation — Weapon  Delivery 
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Figure  14b.  Standard  Deviation  of  Tracer  Error— Elevation— Weapon  Delivery 
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Figure  16b.  Standard  Deviation  of  Tracer  Error — Azimuth — Weapon  Delivery 


MERN ( RflD ) 

-0.04  „p.00  0.04  0.08 

_ I _ tn_i _  i  i 


MERN( RflO  ) 


ELEVATN  TRACER  E 
SU8JECT  33 
TRAJECTORY:  ZIGZ 
CASE  40132 

- EMPIRICAL 

....MODEL  PREOI 


CO 

10 


De\ 


Figure  18b.  Standard 


AZIMUTH  TRACER  ERROR 
SUBJECT  33 
TRAJECTORY:  ZIGZAG 
CASE  40195 

- EMPIRICAL 

....MODEL  PREDICTION 


AZIMUTH  TRACER  ERROR 
SUBJECT  33 
TRAJECTORY:  ZIGZAG 
CASE  40195 

- EMPIRICAL 

_ MOOEL  PREDICTION 


ELEVEN  tracer  error 

SUBJECT  33 
TRAJECTORY:  JINK 


CASE 


40132 


— empirical 
_ model  prediction 


V 


Figure  22a. 
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Figure  24b.  Standard  Deviation  of  Tracer  Error — Azimuth — Jink 


Section  VI 
CONCLUSION 


This  report  summarizes  the  modeling  of  a  gunner's  performance  in  a  complex 
AM  tracking  and  firing  task.  In  this  task,  the  gunner  observes  not  only 
the  tracking  error,  but  also  the  miss  distance  of  tracers  from  the  target. 

A  gunner  model  is  proposed  here  which  consists  of  a  reduced-order  observer, 
a  linear  feedback  controller,  and  a  remnant  element.  The  Average  Approxi¬ 
mation  Method  is  used  to  solve  the  closed-loop  delay  differential  equations. 
The  parameters  of  the  model  are  determined  via  a  least-squares  minimization 
algorithm.  Computer  simulation  results  show  that  the  model  predictions  of 
mean  tracking  and  tracer  errors  are  in  close  agreement  with  empirical  data 
for  several  flyby  and  maneuvering  trajectories.  These  results  demonstrate 
that  the  proposed  model  is  an  accurate  and  efficient  model  for  representing 
the  gunner's  performance  characteristics  in  an  AM  weapon  system  in  this 
tracking  and  firing  operational  mode. 

This  human  operator  gunner  model  has  been  incorporated  into  the  MTQ  series 
of  POOl/OBS  AM  engagement  models  and  is  designated  as  program  POOl/OBS  3/6. 
It  is  this  final,  composite  program  which  is  intended  for  use  by  aircraft 
survivability  in  performing  their  weapons  effectiveness  studies.  Documen¬ 
tation  of  POOl/OBS  3/6  is  in  preparation  at  the  Air  Force  Aerospace  Medical 
Research  Laboratory  and  will  be  distributed  separately. 
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APPENDIX  A 

PROJECTILE  BALLISTICS  AND  INITIAL  CONDITIONS 


PROJECTILE  BALLISTICS 


The  projectile  ballistics  equations*  for  an  AAA  weapon  system  used  in  this 
study  can  be  described  by 


R(t)  = 


v  •  t 
o 


1  +  t  •  K1  +  K.2t  -  /K^t  -  K^t  +  KJ 


ot(t)  =  sin  1 


cos  01T,(o)  sin  a  (t) 

_ IB  _ o  _ 

r  -.1 

2  *  1  2 
1  +  sin  a0(t)  ~  2  sin  a^(t)  sin  9]_g(°) 


aQ(t)  =  0.001  (K6t+K?t  ) 


with  constants 


v  =  930.0  m/sec 
o 

Kl  =  0.19036681/sec 

K2  =  0.14851772/sec" 

K3  =  0.02565365/ sec^ 

K.  =  0. 02274266/sec'; 
4 

K5  =  0.01660118/sec" 

K.  =  5.202815 
6 

K?  =  0.4863915 


*See  Milenski,  J.,  "Methodology  to  Characterize  23  mm  Projectile  Trajectories 
for  Monte  Carlo  Simulations,"  Technical  Memorandum  to  J.  Bode  et  al., 
Braddock  Dunn  MacDonald,  November  24,  1975. 


a(t)  otQ (t)  cos  01B(o) 


=  O.OOl  (5.2C2815  t  +  0.4863915  t  )  cos  6  (o)  (A5) 


On  the  other  hand,  within  the  range  of  interest,  it  can  be  shown  that 


K2C  "  A  /K3t  "  K4t  +  K5  -  0,028  <<:  Ki 


We  car.  approximate  Equation  (Al)  by 


R(t)  « 


v  •  t 
o _ 

1  +  Kjt 


Thus,  for  a  given  range  R,  we  can  solve  for  the  time  of  flight  from 
Equation  (A7) 

=  R _ 

T  ~  930  -  0.19R 

INITIAL  CONDITIONS 


In  the  computer  simulation  program,  the  gunner  starts  firing  at  t  =  tQ  -  T 
with  the  first  tracer  round  reaching  target's  slant  range  at  t  -  z  .  The 
gunner  model  is,  therefore,  turned  on  at  t  =  t  .  The  initial  conditions 
used  in  the  gunner  model  are  computed  as  follows. 
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Azimuth: 


02B(to_T)  =  "°-025  X  sgn  [®2T(to“T)] 

82B^to)  =  _T  62T(to') 

W1  (to)  =  [02T(to)  "  02B(to)]  C0S  8lB(to) 

W2(to}  =  [02T(to)  "  92B(to-T)]  COS  0lB(to) 
w3(to)  =  W4(to)  =  ° 

Pll(to)  =  P22(to}  =  °-0002  X  cos  0iB(to)’  P33(to)  =  P44(to)  =  0 

Pij(to)  =°  V  i  *  j,  i,j  =  1,2, 3,4 


A 

where  T  is  the  anticipated  time  of  flight  of  the  tracer  fired  at  t  -  t  . 

o 

This  value  is  normally  available  in  the  engagement  model  as  the  difference 
between  the  time  of  intercept  and  the  time  of  fire.  In  general,  a  rough 
approximation  can  be  obtained  by  using  Equation  (A8)  with  R  =  R(tQ) . 


APPENDIX  B 

A  HILL-CLIMBING  OPTIMIZATION  ALGORITHM 


Let  a  be  a  vector  of  dimension  p  *  1.  The  cost  function  J  in  Equation  (30) 
can  be  rewritten  as 


min  J  O) 
a 


=  mm 


a  j  - 


Z  /  {[c"S(t)-i'j 


(t) 


+  £ 


-1 


(c  s (t)  -  Sj  (t) 


dt 


(Bl) 


The  iteration  procedures  which  adjust  the  values  of  a  in  search  for  a 
minimum  of  J(a)  are  described  below. 


1. 


Choose  an  initial  guess  a^  of  a,  and  e  which  is  a  termination 


lower  bound.  If  the  absolute  sum  of  changes  in  a  within  one 
iteration  is  less  than  or  equal  to  t,  the  iteration  process  will 
be  terminated. 


Compute  J(a  )  and  call  it  J  . 

— o  o 


3.  (a)  Choose  a  set  of  p  orthornormal  vectors  v_. .  cumulative  vector 


f. . ,  increment  coefficients  e.,  cumulative  increment 

-zi  l 


coefficients  d.,,  success  indices  q^,  i  =  1,2,  ...»  p. 


Initially,  let 


v.j  =  =  (1,0,  . . .,  0) 


v2  =  i2  =  (0,1,0,  .  V 
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Vp  =  ^  =  (0,  .... 


T 

0,1) 


e.  =  0.1 
i 


d. 

i 


=  0 


qi 


9 


i  —  1,2,  •*.,  p 

(b)  Set  parameter  axis  loop  index  i  =  1 

4.  Adjust  the  parameter  vector  according  to 


a  =  a  ,  .  +  e.  v . 
—new  —old  l  — l 


5.  Check  constraints  on  a,  e.g.,  if  £(i)  are  constrained  to  be 
nonnegativ  then  set 


a 

—new 


(i) 


a 

—new 


if  a 
—new 


(i) 


<  0 


6.  Compute  J(a  )  and  call  it  J.. .  Then,  do  the  comparison  test 
—new  1 

If  -J,  >  -j  ,  set  J  =  J,,  proceed  to  (7);  otherwise, 

1  9  0  1 

proceed  to  (9) . 


7.  (rj  Store  accumulated  increment  along  each  axis 


d  .  =  d  ,  +  e . 
xii 


(b) 


Set  a  ,,  =  a 
—old  —new 
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»*« 


mPSBIQPtY^Hh 


Rotate  orthormal  axes  v±  according  to  Gram-Schmidt  orthogonali- 
zation  procedure.  First,  store  cumulative  vecor  1. 

+  d2v2  +  ...  +<1^ 

i2=  d2i2+"-  +  dP^P 


Second,  compute  a  set  of  orthogonal  changes  3.,  and  new  axes 

V. 

— 1 

®1  =  ^-1 

ix-Trtrr 

*2‘h-  {k* 

h 

1-2'TT^nT 


MMMflMMM-liMMM'IMNHMt'MrtMPhMMM  (MMMMMMHM  FHNKM  MMNM  MNNMMMrIMM  MMMMMM'IMMMMMMMMMMMMMNMN'I 


W6W,T  ?5G,CH79-3(H.-»— L-7-frG255,WEIrgSa-3960 - - 

MAP ( ON ) . 

GGMME-NT-, - »  -  PI  Hfc-VARY-If^.r-J-.-Q-.F-. - 

COMMENT  .  »•  ElFIT,  I0=L76Q295,CY=3?  TPETPD0T=ALPHA(1)*THET8(T-TAU>*» 

FT  N. - 

ATTACH,TlPcl, 3AT  M6RAZ,  CY=3, MR=  t • 

ATTACH,  TAFE2',-a4TM&4tAG€fi~rW^Hfis4T - 

LOO. 

- P kOG k  A-H— OPT-(-TA P E-lr T  A P E  3 ,  i  NPUT-,  C-U7PUT) - 

DUENSION  ALFhA  4  8),  PSI  (8)  ,A(6)  ,0(6)  ,  ALPHAN<3> 

- — 01  ME NS-T-3N — STrF5~t"7")  ,£  (8 4  6,8),  Z 8~,8 ),  F-T7 ■)-, F7T 7-,-8) — - 

CCiMON/AkRA  Y/X  (1LOO),S(1QOO),  ELOD  (1000)  ,TX(10  jQ)  ,TS(  103  3),'?  AN  (1000 
- -t) - 

gghmon/5/tau,cq,oel,nstp,ndim,to,ipt,yio 

- bOGIG-AL'-PHr^F - 

INTEGER  Q,SUEIT,FR 

- -LG=0 - 

00  1  1=1,8 

- 00-1-0  =4-,-9 - 

V  (I,J)=<I  (X,  Jjs-G  .0 

- iP.  (I  .SVJ,  -V(i,  i»*XHtr-i»-trO - 

Z(I,  J?=i.O 

- 6-<I-HWr6 - 

1  C  CNTINUE 

- FR»1 - * - 

READ  *,  Kl,  NDI M, T  u  ,IP T»  Y10 

- PRINT — ,arKl,Ki--I  M  rT-g-ri  ^-Tt  Y10 - 

42  FC^MAT (1H1, "AC. OF  PTS  =  ”,  In , 2X  , "ORDER*  “,X2, 

C.5) 

- KT  *T-t3/-thr8-3 - - 

K=K1-KT 

- F?6AG(-1:t-43)  (OLh6  »  AZi-AZO,  AZGi>TfeL  t€LO,aOHl-,T,A7MN,OUM2rA?SOfGOM»f - 

C  1=1,  KT) 

- tjO  ty4 — I-^tTrK -  -  - — 

REAC(l,t3)  CLM6  j  AZ, AZD, AZCO,  EL ,ELO, 0UM1 ,T , AZMN, 0 UM2, AZ5 0,CUM3 
- i-f-<  EO)Hl9-)-45r45 - -  •  - — 

43  FORMAT (12011. 5) 

45 - I-r-T  *-«■♦,£  ,0)  GO  -TO  44 - 

1H=(I-1)7IPT*1 

- g  LOO  MMF^-GU-Ml - - - 

X  (IH)  =  0'JM2 

- S(Ilt)  *  ODM3 - 

RAN(IH)=7 .5 

- lF(OOMO^tT',  ge77'r)  RAN  ( III)  »OUHC./(93C'  i-  >19»0DMC) - 

44  CONTINUE 

- 0a-47-T~»lrK" — — - -  - - - - 


* 

2 


i 


* 


i 


k£AO(  2,  -,6)  J  IF£,T  AZ,0UM4»7AZSD»  CUM5 
- IF(EOF-<i)  )Hfc,-»yJ - - - 

46  FORf-.AT  (5G12  .5) 

40  - H‘*Fr4FrHr4»£T^  ^  TP,,ft7  —  ■■■  — 

iG=(i-u/iP  :*i 

- TX(  IG)-»-JUH4 - - — - 

IS(IG) =  3UM5 

47  - CONTINUE - — - 

0  £L=  J . 0  3  * IPT 

- b~0«  1  t-5  4 - - - — — - — — - 

N£TP=IH-1 

- k£4£r*  H  4i;PH-/  r8) - 

R£AC  ‘ , EJS 

- A  J=C-t€ - 

1L  00  8  1=1,8 

£(!)*. 1 
0  (I)=C  .  3 

- A1IM2  ,  - 

8  CONTINUE 

- m=suc-it=-: - - - 

oALL  INTO  (ALPHA,  AJ) 

- GLOU — A-J - - - - 

PRINT  -43,  AL  (HA 

4U.  — F  CRMAM"OALPHA=— ^TC.12^5rS-«-VLrC-l-5-.&rIi^lit^-lSr&r‘1PiV&14»-5r:l^-» — 

C  G12.3,''i",G12.3»“l“,ol2.5,”t“,G12.5) 

- PRI — ■*  1 ,  A  J  .-I  T-,SO(~IT — - . - 

41  F  CRMA1  (  "  J=  ”,  G 1  2  •  5  ,  ”  ITERA  flCNS*  ”,  1 5,  “SUbI  T  £P.A  TIONS=  ”,I5) 

■11— 00 — T3L -if-8 — - — — — ■ —  — .  — ■ .. 

SU5  IT  =  S  J[  IT  <1 

- 00-12-  L=l, a - - - 

A  IPnAN  {•_)  =  ALPHA  (  L)  «■£(<)*  V  (K,  U 

—  1 F  (  ( u  • ’j-;  • -3 )  *  AN  0  .1 L  •  Lt.5)  3-  GG—T  G-12 — - 1 - - - - 

J.F  (ALPHA  tl(L).LT  <C  •  j)  ALPHA,- HLJ--.-ALPHAN  (L) 

-It; - CONTINUE - — - - 

CALL  I.NT6  (ALPhAN,  *.J) 

- IF — ( AU-»-jT-rOL£-J  3 — GG — TO — 20 - - 

CLDJ=  A  J 

- -  PRINT— *«3rAL-fhAN - - 

PRINT  41,AJ, IT.SUBIT 

- - - - - - 

£  !<>  =3  ‘£  ( K) 

- iiO— 1-3— M=T,  8 — — . . . . . . — — — — . - 

nLPHAMMALFHAN(M) 

15 - CONTINUE - 

IF  (A(K).oT  .1.5)  A(K»  =  1.0 

- G-C-TO-2  j - 

2G  £(<)■=-. 3*  £ (  K) 

.  IF-(A(<>. Lc-rtr^-) — A4-K>-=4.-P - 

25  C K3  C  .  G 

- 00 — 33 — L-»1t3— - * - - - 

xF  (A(L).LT ,C.5)G0  TO  30 

- GK=4-r£ - - 


I 


t 


( 

: 

■? 


88 


J 0  CONTINUE 

- I-f— (-€K-r+)E-»  0-»  C)  &-6 — TO— l-M - 

SUM1=SUN2= a  .c 

- DO  -32-Mat-r6 - 

suNisSUii+xni,*)*^ 

- I  <5,  - - 

32  CONTINUE 

- X  lsS-aR  HStiMH - 

X  2=S3FT  (SUM  2) 

- x-:*=  X  tv-x-2 - 

PRINT  33,OLD  j,Ai_PHA,Xl,X3 

33  - HHHA-T — (u  tU»  - - rO-15  t  3/'  Ai_P‘ii*-»t~i  G1  £?<  Or  » *T&t-?-T^T*<-»T^t-gT?r 

2v**»612.5**li,i612.5*‘:»*»Gl2.5t->“»612.0»',l9»Gl2.5»*)*/ 

- 3«~— X  I-Hr)-=  ■■  ,  l-12  i-Olg-rSH - * - 

GC  TO  110 

10  0 - GGNT-I-NUE - - — — - 

105  GO  TO  11 

110 - SG43«l-»-) - - - — - - - - — — . . 

DO  115  Lsl, 7 

- OO-tl-5— Isi-r-t - - - 

F  (LMC  . ) 

- *2<L,M) -O  r# - 

115  CONTINUE 

- GG-l-t-7— J»-lrr6- - — - * - 

CO  117  <s 1 ,  i 

11 7 - XlN  JrK)-C-r0 - 

DO  12.  1=1,8 

120 - SUN  3  =  SUN3+4-ieS<^-(-I4-)-) - 

IF  (3UU.LE.EFS)  GO  TO  1000 
CO  130  N=l, f 
DO  13k  Jsl» c 

- ^Nr)|sa(NHf«MT - 

13G  CONTINUE 

- DO- 1-4  f— J  - 1-,-6 - 

CO  14C-  L=  1, 

— - DO — J,  (  —  ■  ■ —  ”  ■  1  ■ 

XI(J,L)=XI(  .,L)*Z«,L) 

140 - CCNIINUt - 

SUM4=  C  .  ) 

- CO -1-51 — J~1t5 - 

150  sON^-b'JNL+XI  (1,  J)**2 

- —  — SUN^-s-SORT-  ('b~U4) - — . — ■ — — —  ■  —  — — 

DO  155  J=l,  6 

155 V-(l-r-J)— XI-t-1  rJ)V-SUK4- - 

K  CUNT  = 2 

16  9 - Q=<GUN-T-i - 

CO  171  <= 1, Q 

- M>-16( . L»lr4 - - - - 

F  Q)-F  (Q)  +  XI  { XCUNT,  L )  *  V  (K,  L) 

loO  -  -  CONTINUE - 

DO  170  N=l,8 

- -  .  F2(Q,H)=FtG  7  -*  V7f'  ,-M) — F2  (  0,-t4) - 

170  CONTINUE 


CONTINUE 

00  1-91. — 1»  1 1  8 - 

o  (<0Uf.T,I>  =  >I(KCUNT, I) -F2  (Cl*  I ) 

-  bUMHO)  =0-.-3 - 

00  2H  Ms  It  8 

— 50M^f&-MSUI+V(4f-*-b4X-€RI>«VH->-*— 3 - — — 

bUM5  ( G )  =5(K  i  <bU.N5(Q)  I 

V  «CUN7,M>  =  (KOUNT  , ID  /SUM5TCI 

KCUNT*.  <0UN7»1 

— HHX<WU-M>-  GG~T<9~1»» - - - - - 

i  i=mi 

SU3IT-0 - - 

F.K=1 

— OC— 23< — <«1H> - 

£.{<)  =  .  t 

— 04Kt-»t.— »■■) -  - - - - — — 

A<<>=2.3 

CONTINUE - 

OC  TO  11 

CALL— tXIT - 

ENO 

— &itd POUfINC — INnj  i -AtP-H-Ay rr»  — .  . —  - 

COMMON/ AF.P.A  i/X£MHluJj)  ,S£H?  { 1  jj.),  £ 0 D  I  1C  30) , TXT  1113) , 

- -lrTAUH-WC) - 

COMMON/ 3/ iT  AU.CL ,t)EL,NSTP,Nu'f TO  ,IPT,Y10 
—  0 1M6NS ION  W<-C>',^<^r,^,ckH-Cr->T-P2TCv^)»ALPMA{»)tA<C,4>- 
l,e  (lo)  , £[ <!£),£.  INT  U6),E  AU,  4 > ,  E AI NT (4 , 4)  , “  {,), C0C(4, 


TS(IPJQ) 


A  >  ,  C  C  +> 


EQ'JIv/A  LENCE  (A(i,l),3(l))  t  (CAdtDtEbtl  ))  t(£AINT  (1,1)  f'3INTd)) 
- OATA-KT/lr/ - 

3~I N I T I AL IZ  A-T-IO  N - 

■> 

■* 

- - 

CO  1  1  =  1,  NO 

- GO— 1— U-lfNQ - 

P(I,J) = J. 

- ^OC(IrU)-0-» - 

L  A(I,J)=). 

NC1=ND-1 
N  C2  =  N0-2 

- NAA-NC-/Z— 1 - 

CO  11  1=1, NO 

- Hflf={-. - 

0(I)=C  . 

- F-«4-=4-, - 

LI  P<I,I>=}.02  0 

. — F-<  1,1)-  IrOO— 22-V-i Z - 

P(2,2)=).JC  0E32  46 

- X  o  ( 1 )  —  *3  ,0  •jV't’ZG - - - 

A4( 1) =0. 

- fcGH4 1  34-1-) - 

£00H<1) =C . 

_  - C  3-0, - 

ARj  =  -C£urALPHA(  2) 

- IF (A3G.GT  .-2-0,^)— S1-Z-XP4AT1W - - - 


-J1"  AmJA  — I.  ■  / 


fr 


c 

-C -GOMPU-TE-ANO  5T0  E  S'!  ATE-S  Xg-AMO  X4 - — — - - 

C 

c  ■  —PRINT — 9? - 

9  7  FORMAT  UX ,  "TIME",  4X, "TARGET  V£L",2X,"EST  VE.  ERROR",  2X,  "EST  TAR  v/ E 

- l-L-'V) - 

r=TG 

- ^4-;-^— l-rKST-P - - 

T=T  +DE  L 

C - if  {-(MOO  T-rOR-^-LKK  *£0  rl->  T-P«I-MT-9>,4 ,  X3<  <L X-4-LKKJ - 

C  1  ,X3KX>-X4  ikK) 

9o - F  CR  M  A  T-  <  ^S-LR-r-S) - 

K1=<K  +  1 

- K2--{  K-4 - 

A3 { Kl> =<3<KK)+£LC(KK)»0EL 

- 1F<  ALFHA<2»-.£G»-v-»)~GG--T-G~4 - 

X M  K1 )  =  si  X4(*Ki  ♦  E'JOCKK)*  (1.- SI)/ ALPHA (2) 

- MiO  TO-  J - - - 

4  X  *♦  C  K 1 )  =  X4(KK)*-EGC(KiO“  QEL 

-3 - GOM+-INU* - 

C 

C-- COMPU-TE— ANO-STO-^-E-EEO-INA-T-ES— T-AP.E-&-T-- VELOOI TY— AMO-AGCF  LERA-T-ION - 

C 

- tGH  {  KK  )  -X  3  ( *K»*-*-M  XK-) - 

IF(K2.3E.l)  EOOli  ( K<  )  =  (  EOH  (<i<>  -£0HU2)J/  CEL 

-lu - G04T-I  NO* - 

NClM=NO 

- N1—N3---2 - 

C 

C  START  -I  NT  EGRAT-ION-fcOOF - 

C 

- - - - 

T=TC 

. H  ( 1 )  =  V  -10 - 

H  (2)  =C  .)2>t736 

-4.MEAN-  <  {(  1)  -XENP<4T5>  M(-2  )  -T-X-(-l-) •>-*-*••  2 — - - 

SS3=(S0TT  (F  ( 1,1*  )-SEMPU2  5>)'  -  2MS0RT<P  <2, 2) )  -TS  (1>)  **2 

0  0—1  9-S—  --l-E-5,-NL-4P - - - 

M=<><+1 

H«.TAO<*:<.-/CeL - - - - - 

RA  =  NAA/7AU'  kk) 

- - A2»1«~HA0<  *I>-T  A1PK+CHV-OR - 

COR=CC4.Tl.p!iA  (3i 

- ^  {  1-y  -  '  — - —  1  -  *'  — — — — 

A(l,2)  =-C.>  ALPHA  (4) 

A(2,NG1>=;AU,1>--AL-*HA<-113-A2 - 

A  12,  NO)  -  A  ( 1 ,  2i '  4  LPH  A(  1)  *  A2 

- GO  2--1—  1tN02 - 

01=1*2 
A ( Jl* I )  --RA 
A  Ul, Jl)=-R/ 

2 - COMT-INUE - 

CALL  LSCRT  (NCIM,fc,DEL,cB,  Ebl NT  »F) 

- GR3-0E-ML-PHAR9-H-ALPHA  M-H-Ag - 

CR4=CC* ALPHA  tb) 


i 

< 


C  (C  (  2»  2)=AlF1-m61~SCAL1 

C  - IF<  <M00<KK,-U---)-»-t&.-3T-rGP-.  TK*-»-t4rlT4-PRlNT-99»T,W-<l-)-fYi-}t  SHF AV) 

C  1  SQfTtP(l,il ),SSC 

99-  FORMA  T^>G  12  tSJ - 

C 

C  -  COHPtHe--HE-A-N--PRA-c*-£  NC-  ERROR - — - 

C 

- 00  -i-1  •.—!-»  1  rN  c  - - 

00  120  J=l, TCIM 

- 0-<I )  -0<-f >-■*■&■  A-(ifO-)A-W-<-9 ) - - - - - 

120  CONTINUE 

1-10 - OONH-NOc - — — - - — _____ 

F(I)-U.-CRh)  X3  (KK)  ♦CR4’1  X4(KK) 

- F(2)=X-3<KKF - - — — - 

IFt«.£T  .  M)  F(2)-F(2)  ♦CR3‘  iXtl  KK-M)  -X3(  KK-M) ) 

- 00-1-3-.— i-NHH - 

00  14C  J=li 

- u-TI-H»fr-(-I-H-E  A-IH-T  <  I  i  J )  ■»  F  ( J ) - 

IsO  CONTINUE 

- rt  <  I  >  =  - 

0  tl ) -c . 

13C - CCNT-INUc- - 

c 

G-COMP-tlfE— ERROR— CrUC— Tt>-HE  AN  RRA-&K-I-NG — ERROR - 

C 

- -S  HE  A  N*  S-IE  A  N  ♦  MP-<  Kl-T  T-'-E  -HH  <  2  )-TX  <  K 1- 1 2  S  >  >-*-*-2 - 

C 

C  COMPUTE- eov AR-I-A  NCg-  HA T SIX - - - 

c 


K  «H-fr)-*-ARSRRQOMK<-lH - 


l  «SCAL 


- IFKR-yGf-M)— ttC-<-2T2->=<-ALPHAH>)-fALPHA<?>-A{jS<EDH(-KK-MH  — 

1  ♦  ALfriAU)  K  S(  EOOHUK) )  )»  SCAli 

- C  ALL  -HUuT  <  E  A  IN  TtOOCyNO IMf  NlrRi-Tf}  ) - 

CALL  MULT  (EA,F,  NCIM*Nl-,P2,lO) 

- O— 3  2  F— f  •-  t-f  N  C-F+1 - - - - - — 

00  221,  J=  1,  KIM 

Pil.  J)=P1  a-f  j/*P2  ' - 

220  CCNTINUE 

C  COMPUTE  ERR OF  OUE  TO  STANOAPO  DEVIATION 

C _ _ _ 


^S3= SSOMSOrTiP { 1,  1) ) -SEMPtKl) )*v 2+ ( I  (P<2 , 2) ) -TS( Kl- 1 25) ) »*2 
nP=K°+1 

It  6 - GO-NHMJE - 

EJ=OELMSHEAN+rfT*SSO> 

- REWW - _____ _ _ _ 


E  NO 

-  rUBROOTINc — HGtrT-lEf F  iL>LltH|rtR) 
OUE.NSnN  Ell),f  (1),G(16),H(1) 


I EMP=C  . 

- GO — &- JaIf  t-l-f-t. - 

IEMP=I£1P+E(j)«t  (II) 

5 - i-l=-tl+ 1 - 

KK=  (<—l)  *'(_♦! 

H (KK) - TEMP 
10  6<KK)=T£MP 

- 1  F-<  M*.  £-•)  *4-»-fi-e?-Uft*- - 

00  23  1=1, L 

- 00-2-3— K-I~rL - 

TE’-IP=C  . 

- £J=* - 

OC  15  J=I , LI »L 

- » £.MF-I-£.1F-*£-<-03*-E-(-I-I4- 

15  11=11  +  1. 

- K  K=-(  <  -1  J-M.-+-I - 

20  H  (KK.)  =  TEMP 

- L-2==L— 1 - 

OC  33  1=1, L2 

— - l-3=  I>  1 - 

00  33  J=L3, L 

. — K 1=  <  I -1-)  *  L-+ - 

<2=<J-1)*L+1 


....  -S03P0U  TINE™  £-S0«.  T  ( KQIK,  A,  DEL-r£A  ,£A  INTrNT) - 

OIHENl  I3N  A  ( 1)  ,  L  A  ( 1  ) ,  £ AI  NT  ( 1 )  ,  C  (EF  (  3C ) 

C  - $£»  j  Ert=£XF(,*-G£L>rSAlNT-lNT-:0«At-£^ -O-M-OEL- 

NCIM1=  N3IM+1 

- ^^4=^•.3  - 

NVU  =  NT-1 

- &C£F(  N  n-i. - 

00  1)  1=1, NT M 

- i-i=Nr-i - 

10  CG£F  (IU  =  DEI-'LU£F<II  +  1)/FL0AT(I) 

G - -NT-^UE-T-pc-A^-E-S-AS^-a - 

GALL  OIAG(N (lh,  i AINT , A ,C0 £F il)  ,C0EF(2)) 

-  00-  6-3 — L— Or  AIT-  -  — — —  — — 

CALL  MULT  (A.  EAIr.T,N!3IM,NN,EA,l) 

- 1 F  {  L  •  £  •  J  •  N  ^)-00-T-0-?-3 - 

60  CALL  GI4G(N0II',EAlNT,£A,l.C,C0eF(L>> 

- 70— -00— 6-3 — I-T— 1-,-N  Kr  M0-IM1  —  ■  —  1,1  ——————————— — 

EA(i;>=£A(II)+l.C 

- 8C— GONT  IftOi - 

ENO 


§u 

Ki 


DIMENSION  A  (  1) ,  (;  ( 1) 

- - 

NN=NJIM -NOI ) 

- UMl^NC-tH— t - 

il  =  l 

- I H  C 1 »  £-1  .-lT-y-)-60-T-6-t3 - 

LC  5  J  =  1| NN«  NLI H 

- K»-J  +NH-1 - 

(jO  a  I=J,K 

- <♦  -  A  <1 )  -  C  l-*6  (-I-) - 

« (II)  -  A  d  X )  <C2 

....  5-  1 1= i I +N  )I  HI - 

KtTIKN 

- 10— uO-  - 

K  =  J  +  IM 1 

-  oO  6  I-  J,'K - 

6  A  ( I )  =f.  (I) 

A(I.  )  i  A  ( I-I)-*- €2- - - 

7  ll=t.T*fOIMl 

- k€TU-?<. - 

ENO 

22  C  ttf  3-J  •  o f  2 y  - —  —  •  —  — 

1.  If  2.,  1.  vs:7f  1.6419,  .  921  33,  .11526F-4, . 14 41 5F -4, , 46165 E-3 


46 


47 


lL 


IF(MT)D(I-I. I P t » .KE.3.)  GO  TO  47 

Ic=*<I-t>/  IPT-n - 

fX(IG) =  0UM5 

- Ti><-I-GH*6UN? - 

CONTINUE 

— u£L=  J,  0-3»  IPT - 

CG=1.28 

—  hSTP-I-H-1 - 

READ', ( ALPHA (1) *  I=1,NPAR) 

READ", EPS 

PRINT  ,EPS, (ALPHA(I),I=1,NPAR) 

—  AJ-»C-.{ - 

CO  6  1  =  1, NP  A  R 


- 6  <I)*r~l- 

0 (t ) =t  .  ) 


- A-H-W-r-] - 

S  CONTINUE 

- C-)=SUH-T=-> - 

CALL  INTGIALFhA,  AJ) 

CLOJ-A-J - 

PRINT  4J,  (ALPHA  (MP>  ,MP=1,  NPAfi) 

4b - r  GJJ 

C  G12. 5 ,  ”*  **»  6  12. 5  »#,l  "* G12. 5) 

—  — PRINT — 41 ,  AU-,*II-»-GU6-I-I - 

41  FORMAT  <••  J=  *\  G 1  2  •  5 » “  HERAT  I  OKE= 

11 - OG-  IOC— X— lrNFAft - 

SU3 IT 1  SJo IT  »1 
— - - G  0 — 1-3 — b«-i-rN  P'AFt- 


,15,  "3U3I T ERA  TI0NS=  **» Ir ) 


12 


ttLPHAN  ( A  FhA(l)+t(K)*V(K,L) 

--1F<  <L  «-6S.-2)-.-ANt>.-tfc-»L:€-.  4-)  1—60 — T-6-1-2 - 

IF  (ALPTANCl  ).LT  .u. 0)  ALPHAN(L)=-ALPHAN (L » 

-CONTINUE - 

CALL  I  .N  TO  { A  LPh  A  N  ,  A  J ) 


OLOJ'-AJ 

-  PRINT -h3»  (AL-FMAft-fM°+-rMP»t-rNPAiR-)- 
PRINT  -tl  t  A  J  i  Ii  »  S  UBI  T 

GK)=CH<>*ETK) - 

E  «)--3>E<K» 

— uO- 


15 

2t. 

25 

3  c 


-16-No-l-r  NFAf; - 

AL*jPA('1)  =  ALPFANIM) 

CONTINUE . . 

IF  <A<K).l>7 .1.5)  A(K)“  1.0 
l-C  TO-2-> - 


£(<>--. 5»£(K) 

- A-f  K  PrU  C-.lrG)  A<  *<•)-■ 

C  K=  C  .  C 

-  UG  -SO-L-U-rNPAl - 

IF  (A(U.LT  .0.5  >  GO 


-ErO- 


ro  30 


CK*1.I- - 

CONTINUE 

iU<il-:SUN2=0  .0 

CO  32  H*1 1  NFAK - 

SUN1--5  UNI  *XI  (it  •N);->2 
iUH2sSUN2»XI  - 


96 


Si  CONTINUE 

- X 1  a  SQft-T  <  SGfFH - 

X2=SQR  f (SUM  ' ) 

- Xi=X-l/-X2 - 

PRINT  3  3 1  OLuJt  (ALPHA(.MP)  ,MP-  1,KPAR),X1,  >3 

33 - FORMAT — IM-d* - <~,-&-l'2-r-5/ 1 — At-PTt  A»(-‘-y-&i-2v5r“  •  ^rGl?  <Fr*  •  *rG"££*‘5r 

2*  »*,  712.5, •*'»%  G1 2.5,*  «l  ,61 2. 5,  >(•', G1  2.5  ,  +  )*/ 


GO  TO  110 

ICG - CONTI  Nile - 

105  GO  TO  11 

HG - SUH3=WT - 

00  115  L- It  NFARt 

- BO— 115-  H»' li-  PAR - 

F  (U»-G  .3 

- F-2a,-H>-=W0-: - 

115  CONTINUE 

- 0  0-4- 1-7 — 0— 1-r  NF  AR - 

00  117  <Ui  NPAR 

117  Xi(J,N>=0.0 

00  121  I* It  NPAR 

12  o - SCN3=  S  UMS-M-A£S<-0(  «->-) - 

IF  (3UN1.LE.EFS) GO  TO  1000 

- GO  -13£ — N=4»-PFAR - 

OO  11C  J*i,NFAR 


130  CONTINUE 

- 00  -14C—  -I— lt-NFAR - 

00  1AC  L*l» NFAR 

— - CO— i-iC — — far - 

XI(J,L) =X1(0,L)*Z(K,L) 

-lvO - -CCN'.-I-NUE - 

SUMA*(  .  J 

- 00-1  Or— J-l  ,NFAR- - 

15C  bUN‘---SUMA+XI  (1,  J)>-»2 

- -SU  1a=SQ  tf-tS-lfM - 

00  155  J-l.NPAR 

1 5  5 - V-(-l  rU T-XI-U  .u/'/  SUN't - 

KC'JNT*  2 

159 --  G*<0UNT-1 - 

00  175  <*1,0 

- 00  -160— L-l» NFAR - 

Fn)-FO>  XI ( KOUNT, L)*V  (K,L) 

-1CT - CC-NT-It*Ue - 

00  170  1*1,  NFAR 

- F2(G,M)~=F-* O-f-'V-l  K,NT — *--F2  ( OrH)— 

17 C  CONTINUE 

. F(Q)=fi-rT - 

175  CONTINUE 

- 00-1  q;— ^-1-rN  PAR - 

19C  j  (<C:JNT ,  I )  =  IUGUNT,I)-F2  (Q,  I) 

-  3UN5(0)-0-.G - 

cO  20C  M=  1,  NFAK 

2b  &  -  3  015(0  »-*SUHS  (C-MTrfKOUNF'rMF* — £ — 
bUN5  ( &  I  sSuRl  (SUMS  (“J) ) 

- 0  G  -  21 5—1=  tr-NFAR - 

V  «0'JNT,!1)  =  ;  (f.0UNT,rt)/3UM5(3) 
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o  o  o 


21  &  — uC'JT  tUUs - 

KO'JNT:  KDUNT  +  1 

-  IF  <<OUNT-.Lt.NF)tR>-GO-TO-15-3 - 

n-ITtl 

- s.  - 

fr-  i 

oC  25C-  -K-1,  NFAR - 

£(<)-.  1 

•  0  ( K)  =  C-*-J - 

MlK)=2.J 

2,0 - *>GN  HNU€ - 

GC  TO  11 

UCO  0  ALL- -bXTT - 

£KO 

bU-TC 0UTINS--  i  M Gt-AT.FHA  ) - 

CC9 MON/ ARRAY  /X£i*F(iaOG  >  ,SEMP  ( 1  3  It)  ,£DDC  13  00)  ,  £LG  (n)  0), 

- l-y-T-X-M-i-}  3  )fT44tannr) . . . .  . .  — •  — —  ■■ 

CCMMON/S/STAl.C-  »CEu , NSTP *NJ » TO  ,IPT,Y10 

. OH£NS  I-Olt  ■  H  1A)-, -PR-*,  m)-,  Fl  H«, «,  )  »P2{‘*t  ••)•»» L°r(A(  6  )  4) - 

i,L(16)  ,£&<16),E.:  iNT  (16)  ,EAU,u)  ,EAilO  (t*  ,4 ) ,  F  ( -4)  ,  COC  ( 4, 4 

-  --  -  2  »-X3  <  lv- 3  i)-,  X-Ml-G-t-  6>,-G9wR  l&Ou)-,  F.OOhlluOS) - 

£C‘JWAu£NC£  (At  1,1)  .3(1) )  ,  <£A  <  l,l)t£0(l  )),  (2AINT  (1,1)  ,£ 
- GATA— fc-T-Z-l-,/ - - — 


TAuunai 


>,:  (a) 
8I»T(1>) 


— - li  G-1-- 1 — 1 ,  NO - 

00  1  J=t,ND 
P<I, 

CQC(I, J)=0. 

1 - A  ( I  ,  J  }  =-3  r - 

N£IV=N$TP/3 

- — N vl* NO— l - 

N  C2= NO-2 

- N*AM3C-/->-l - 

00  11  1=1, NO 

- WID-'C-. - 

C (I ) =0  . 

— F  <1 )  »(i  v - 

ii  P(i,i)=).:3  c 

- (HI  r  lr)= 9-rrt — 1  <*9-3-? - 

X  i(  1)  =  -  J  •  kG  202 

EIH(l)  =<3il 

- - EC-3H<  l)a-3-. - 

N  £  I  M = N  Q 

- Nl-»NO-« - 

Cb=C0S (£LG( 12b) » 

W(2)*-3, £13351^06 - - 

H  (1)  =Y1  J*  08 

P(lt  l)a-»(l,  1>CG  -'-2 - 

F  l£,2)  =  ( C  . 0  JCtE-V  C8)“'  2 

bS3=(iO?T(Pli,U  )/C8-SEMP(lJ)~*2MS3RTCP< 


2,2) >/C3-TS(l> )*«2 
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- S  1*0  T - 

00  13  K<=1,  NSI P 

- -M=«*t - - 

CB=C0S (£LG(  •  K) ) 

- ARG »— 5 £•.»■»  LfH-A  ( 1  )*C0  —  -  —  — -  ■  ■  -  — —  —  .... — — —  -i  — 

IFCARL.GT.-2C;.>  S1=EXP{ A  F6) 

- X 3 { K 1)  »  <3  fK-K  )+£V&( KK-F*-B£L - 

IF  <  ALT  HA(  1)  .EC. I.)  GO  TO  4 

- X4 '  K 1 )  »:> l^-f*4<< X-)-*F. 3d+KIO-a8  EL-) - 

GO  TO  3 

4 - A-<*  (  K-lH-«*4-f*V)  ♦E-ffriKK) «  OCL - - 

3  CONTINUE 

- E-OH  <  <il-X-3-t-K-13  -y-4  <  K"H - 

IF(KK.&£.2>  EGDt-.  (KK )= (  EOH <K<)  -£QH (KK~i)  )/DEL 

lu - CONTINUE - 

I=TC 

- oO-±jt—  <K*l-g  it  N5TF - — 

Kl=<<+  1 

- M=T  AU<  K-Ov-8-Ei.- - 

RA=NAA/TAU( KK) 

-  00  2  I*l»-NB2 - 

-1=1*2 

- A  (-J  1-,-I-H^A - 

A(J1,J1)“-R  i 

2 - CONTINUE - 

TFSOn= (ELG(Xl)-ELG(KK) J/OEL 

- Cb=COS  <€LG<  kK)  ) - 

Tfes-ThElO'T  /MEL  G(KO)  ! 

- sCAb-H-GT^-GU-)'**' 2/ GEO— . . . — — - -  - - — — - 

T=T*D£L 

- A  2»  1 .  -H  F A  U  (  K  i,  -  T  AO  (  «K- W  ✓(}  Et - 

CuR=CC*ALPhA  (21  C6 

-  — - A  ( 1 ,  1 1*  -C-Oti  + 1  - - 

A  (It  2}  =  -Ci»  *  ALPHA  (3)-xCB 

- A  (  2  *  nC-H-»-0".-fv  A-E - 

A  (2  ,  NL)  =A  (1 , 2)  <  4  2 

- A(2f2)-»T-6 - 

C 

C  -  COMPUTt-FF  A  N5-FF  3-C-N-HAT5U  X--E  A  -  A  N 9— I-T-S-I-N T£~G,VAL— &«-NF - 

C 

CALL  OSCRT(NCIM, fc.OEl  ,Eu, C0INT,5> 

CR4=CC  *4LPH A  (4) *  C6 

C - 

C  ST ART  INFEGRATIGt'  LOOF 

C - 

C  1F((M00(KK,  .EO.Gi.OR.(KK.EG.l))  PRINT  99, T, W( 1) ,SM£AN, 

C- - L-SQRT-(-Pfi-t-ibH  - 

99  FORMAT (5312.5) 

C - 

C  C0MPU1E  'IE AN  TRACKING  ERROR 

C - 

00  110  1=1, NCI M 

- 0O-1-3G  i«lrN  SH» - 

u  (I )  =0  1 1)  *E  (I,J)«H(J) 


99 


120  -CONTINUE - 

110  CONTINUE 

- KU)>( C-3-C^ ) -  X-3-T KK-T-KTRC*  XVLKK-T - 

F (2) =X  3 (KK)  CL 

- -i-F  («-,0-F-.-M->— R  Kli  M>  X3  (  KK  H))JA2 - 

'JO  13c  1=1,  NCIM 

— —DO  !<♦( — -J3-!,  t - 

G(I)=G(I)  *EA1M  l  ItJ)*F(J) 

lvO  -CONTINUE - 

iHI)=G(I) 

- U  <1  )=.R, - 

130  CONTINUE 

C  COHPLTE  ERROR  OUE  10  MEAN  TRACKING  ERROR 

SUSAN*  SNEAN  ♦  CM  1  l/CU-XEMP  (K1 )  >  2*  ( w  {  2)  /C!?-TX  ( K1  -125 ))  *’  2 

- - — - - 

c  compute  covap.ianCc.  matrix 

c - 

CCC(1,  1)  =  (  ALPHA  i  5)  ♦ALPHA  <  6)  A  l-  S  (EQH(KK)  )  ♦AL3HA  (7  )  *  A3S  (EDDHKO)) 

- i-*SGAl - 

cGC(2»2)  =  ALPHM(;)'*SCAL*A2**2 

- i-fRXK  »GT-rHf--C-ftC-r-g-r?-)-M  AL-R-A-fG  >AAfc-PHRP6->— A-GSf-E-DH(-KK--H)-)-*-A-LPH)M-7-)-!> — 

1  /35(tOJH(KK“M)  J  )"5CAL'*A2*V2 

•  -  -■ - CALL — NUuT-1  E-fr  I-Ni-,-CGv  rN-0TH,-Nl-rP~l'T~l-0  )- . . .  ~  —  . 

CALL  N'JuT  (EA,P,  NCIM, N1  ,P2. 1C  ) 

- uO  -22C — I  ,-NGirM - — - * - 

00  22C  J=l,  NCIM 

- f-<  IrFl  1-H-rCM-P-a  ( Iri) - 

220  CONTINUE 

c - 

C  COMPUTE  ERROR  DUE  TO  STANDARD  DEVIATION 

SS0=S$JMSQ  -T(P(  l,l)>/CG-SErtP<Kl>>**  2  «•  <  SORT  <  P  <  2,  2)  >/C3-TS(Kl-125>  ) 

- A«  - 

li>  0  CONTINUE 

- EO=0'L-*-<SME~M-H-r-SS-3J - 

RETURN 

- - END - 

SUBROUTINE  FLLT (  E ,  F ,  L »  L 1 » H  ,  MR ) 

- U-IME-NSIUN— 5-  T-l-F-r  P-Ll-PrC-dO)-,  HM-) - 

DO  10  M,L 

- 11=1- - - 

00  10  K  -1 1  L 

- T  fcHF=C , - 

00  5  J=I,L1,L 

- I-£MP=T£HRA-g-<-dR-P-<-I  I) - -— 

5  11*II«-1 

. -KK=<K-l)M.-*-I - 

H«K)  =  TEMP 

-10 - MKK)-T€M» - 

1 F ( HR . E  J . 1) R ETUR N 
uC  20  1=1, L 
DO  23  K  :  I »  L 

- T€.NP=Cr, - 

II=K 

- fc  C-l-5-  J=I,  L  1-,  L  — —  — — — — — - —  — - 

TEMPs  T£MP*G<J>'  E  (II) 

1  & - I-I—I-I+L - 

KK=  CK”1)*L*-I 


! 


NT  M'JST  G£  AT  LEAST  3 

■  ~CALL~tI-4G  <N-_-iMrtAI  NT-,-A1C&£F+i-T-rGG£f^2-H- 
00  60  L-3  *  NT 

IF (L . EQ.NT) 6C  TC  73  ^ 

63  -CALL  - 0 1 AGT  NO-IKr-t.-AT-NT-»-E-A»-l-»5-r G0EF-<4)-4 - 

70  0C  8 j  11=  1,NN,N'JIM1 

-  fcA(lIT='A«-I4^1rrO - 

oO  CONTINUE 

- 6  NO - 


SU3RJUTINE  LIAG  (  NCIM*  A  r£  ,  Cl*  C  2) 
--  -  JlMENSI}N.-A4-U-,-t^-t> - 


N  lI Ml  =  N  31 M  +  i 
- NN=N-jIM-NOI->- 


APPENDIX  E 

LISTING  OF  COMPUTER  SIMULATION  PROGRAM  FOR 
AN  AAA  TRACKING  AND  FIRING  TASK 

Ho  H  ,  T  2u  iCM7->-J)0-. — L7cv  2-<s5tME-I;,-25  b  30  6fa - 

COMMENT.  -  SlMUto,  IC  =  L7fc3295,CY*lJ  AA4  M00Z6  SIMULATION  PROGRAM  *** 

A I  T  AGH  i-T  Af  cl-r-TA-fMfrR^  Z  ,  1 0  »L  7  b  "23fa-,-0y----3r*R«t-t-« - 

FT  N. 

LOO. - 

PROGRAM  ELF  : Ti I NPLT, OUTPUT,  i AP£l) 

-  COMMON/ i/CG  «  2)  ,u£L,TH,NJIP.  VU-  (E),X  J(  2>-,£L» ELOO,  AZOO.HTAU.RA, A2,NO— 

A1.N02, N)A ,UEL, UA Z, ELTR , AZ TR , I S ET, Zl, Z  2 

G— - - — 

C  THE  PURFOSZ  OF  THIS  FFCGRAM  IS  TO  SIUULATE  AN  E.EVAT  ION  TRACKI  (G 

C  TASK-  )N  -THE— T-RAGEf  -uIREGTEO -F  IRE—  ( *Ciic  -*>)  -  SYSTEM - 

C  INPUTI  THE  ELEVAIIO*  (ZL)  L  AZIMUTH  (AZ)  ANGULAR  ACCELERATION  )F 

C  TARGET - - - 

C  OUTPUT*  MEAN  ANl  STANC  OEV  OF  lAG  ANGLE 

C-**^-*— AGL— AS6LE ■a-* :~:c— I-N-UHI-i-S -OF— AA St AH  -**•  >— * - 

C  TAUI  CELAY  IN  SECONDS 

C  ALPHh  I  PAR  AME  TER -VECTOR - 

C  ELC  R«l  MEAN  EL  LAG  ANGLE  (I.E.  TARGET  A NGLE -PARREL  ANGLE) 

C  AZERkI  MEAH-AZ--1  At- -ANGLE - 

C  ELSOi  STANDARD  LEVlMION  Or  ELEvATICN  LAG  ANGLE 

C- A  ZSQ«— ST-ANtHRO— fa Evt-A-HON-OF-A  Z-LW-  ANGL-S - 

C  ELTSt  MEAN  EL  TKAllR  ERROR  (TARGET  ANGLE-TRACER  ENDING  ANGLE) 

C  AZTRt- MEAN-^Z- TRAtrR-ERROR - 

c  elbar:  mean  el  At  kel  angle 

C  OEL  «  -  T I  ME— 3  TE  P-USZ  0 — I+f-7H£- 1  NT  E‘o*A  T  ION  ROUT  INE - 

C  Y1G  ( 1)  t  IMTIAu  GUESS  OF  EL  LAG  ANGLE 

C~ Y  1G-2 » -TN  I-7-E-Afc— 6-U  Ef  E  -OF  -AZ-  LAtr— ANGfaE - 

C  UEL  >  EL  CONTROL 

C  UAZI  AZ -CONTROL - 

C  C0(l)i  EL  RATE  CCNTRCL  COEFF 

C  CC  (  2 1  I  AZ-RAtE—C-CNf-xOfc— COEFF - 

C  KH  NC  OF  POINT  i  IN  THE  ENTIRE  TRAJECTORY 

G —  K I — NO  —OF — FOtNT  5 — AF-Hrn — T-HZ — —I  RE  t — 7-F-AGE-A — ROUHG — IS — F I  "C  6 - 

C  ELDS  I  EL  ANGULAR  ACCELERATION  OF  TARGET 

C  AZOO »  -  AZ — AN-jUL-AG— AGGELERAT  ION — Oi — t-ARGErT - - — - - - - - 

C  X3  ( 1  >  t  EL  ANGULAR  VELOCITY  0~  TARGET 

C  X3(2)  l-AZ- ANGULAR  vElOC  IT  Y-  OF — FA-RfaZrf - 

C  X*.l  lSTIMATION  lRNOk  CF  ANGULAR  VELOCITY  OF  TARGET 

C -£L I—  cL— A-tt-UtArfa-  t-Or-r-r-lC-fa — OF — rA-RGE-T -  - - . 

C  H ( 1 ) t  MOOEL  PREDICTED  LAG  ANGLE 

C  H ( 2 )  l  MODEL  PRfcJ K f r C  T RACER-  ERROR - 

C  P(l,l»t  VARIANCE  L-F  PClJICTEO  LAG  AFGLE 

C  P ( 2 *  2 1 1  VARIANCE  uF  ARE  01  CTEO-T^AE-E-fi- ERROR - 

C  TCI  fhE  INITIAL  FIRINl  time 


READ  ' »  Kl» TG.  IP) 

PRINT  }fYi,ie.IP1 - 

FORMAT  (lhl,"NC  l.F  PTS=  • 

C"«  12»  "  POiN  >’V  |  - 

K)  =  T  j/ D.C  ( 

— - - - 

i  =T  C 
llfal  M*  A 

IFRINl  -2C/I «J 

Cfc(  1>*  l»?«  '  . - 

L  0  (  2)  :  1.2  *» 

-  NCI*  Nw  t  M-l  ‘  -  - - 

nc:=nuh-2 


’» I-,  2X  ,'*INIT  TI  ”, Gi2  .S'/IX,  "RE'TO  EVERY 


102 


fVaji  cvti* JbLu>. 


I 


. . NAA-NGI-V2-4 - 

0  TAU=7  •  5 

- - 

Ll.SDs.a3C2*  C.5 
AZSC= «  0  3C  2*  C.5 
21=0. 

- 22-G. - 

I  i£T=C 

- GEL-*-}-. - 

UAZ=0 , 

- PRE-NT — ? - 

7  FORMATJIK  ,  2X,"1  IMEM,aXt"EL  V£  L",  9X,  “  EL  ERR  *V:  <»  "ELSD”,  6X, 

- - -i"-£L-GW!1,6XT-,VZ— - 

2| cX j "t L  TR",6Xt~AZ  7R*V) 

-  —  R  £ AD (1 »  2Ft  OUM1  r  A  Z-f-A-ZO-rAZOG  *  c-fc~y£Gt3-r£t.~GGvt  »~AZH>t-t'XvZS&tS  FT- 

C  ) 

- 3G-5- I-s-irrK - 

READ  { 1 , 2)  DU  M.  AZ  »  AZG.  AZOO,EL , ELJ.ELDC.T  ,AZH‘T,  X,  Z  SOtS 

- if  (cOf-l  1-m-r-l - 

1  iFlM0j(I-li  I  PI)  «  NE. 3)  GO  TO  5 

- JH=C-n/-IPT-+l - 

TsTCf(H-i)  CEL 

- T-AUs  T-,5 - 

IF(DUM1.LE.287;  .)  T  AU=0UM 1/ ( 9  3D  .-. 1 9*  OU  Ml ) 

H  1AU=T-«U/GEL— - - 

kA=N4a/ TAU 

-  --  *2=1. -<  TAU-OAUF/GEC- - 

0T4U=T AJ 


- - X-3^  1  )---e.LG-— — ■ ■■■ . *•*  — — —  ...... .  — —  . . 

X3<2)= AZO 

•  -iF((It-U  .GE-.MTAU) -GO-TO-9 - 

TAUR--T  AU 

-  . —I  FT  GUM  1-rL- E  c-AFc-f-FAUR*  ?AO*-AMAXH':-.  fc-,  DUPi/Gfc  3  0-rF- - 

Y!3<1)=-TAUR*ELG-.Gj1''  (5.2*  TAUR  +  .<.85»  TAUR-  *2 )  COSIELO  .  95) 

- ¥-19^*~3ti 0-25-  SIvMl.  ,AZ-3) - 

IFtlH.NE.lJ  GO  TG  10 

- fcL0=EL - 

AZO- AZ 

- fcl9=>  ViGT-H - 

E23-Y1G2 

- b  g--;  o-i-i - 

9  ISET-=I3ET+1 

- ifE  I-a  E-T-.  N  £-.14  -0-3  -  -TG— 1-3 - 

Zl=fcL-<SL :-£i;j  +  .C01' <5.2*TAU*  .486*  TAU*'  2 >* 305 (ELCO.03) 

- <AZx--££C )-  - 

It  CALL  033El6(Ei_£F.R»£LS0) 

- ee^AR— Sr.-Ei.-6PP - 

1FJCIH-U.L  .MTAU)  Y13 (2) *Y1C 2*C0S(EL0AR) 

- GAL-C— GGSA-Zo-FAZ  r  A  ZSG*  EL  BAR  ) - 

IF< <M00<IH-1,IPRIN7)  .EQ.OJ.OR.  (IH.ECUl))  PRINT  b,T,X3tl)  ,f.LERR 

6  FORMAT (llbl2. 5) 

-5 - b-GN-H-fWr - 

2  FORMAT (12G1 1.51 

—  - - STOP - 

END 


\ 


SIHROU  r  INc-  -C-eSEt-fc  <€-LE4Ri  E-LSOF - 

uCMMJN/3/CC-  (  2J  ,  l,EL,KK,NO»  Y  U  <  2  > ,  X  30  ( 2 )  ,  EL  ,EDD,  AZ  OD» M»RA f  A2 

ArNO  1 ,  h  0  irH A-*-,Ur3  AZ  ,  6i-T  4 r  A  £ T  4. r  I  9  ET^-E  1 T l  S - 

DIMENSION  HH.),F(-m),=1U,-)  ,  P2  ( A,  *♦>  ,  A  LPH4  <  8  )  ,  A  <4,  4) 

1 1  d  ( 1  ”' )  (  £6  t-1-64  it:  INT-{164|-€-A<  -«->A)  |E  AI  NT  <*•**♦>  »—(>«>  r  COC  <4  |-*r>-p-0{4> - 

2»  XT  f  It  J>)*X4(1C  >.«.)*  £OH  ( 1 0  C-  C  >  »  EC3H  (luOG) 

£G‘J1  VALENCE— (A<lriH*-<  H-H-<£At  IriHEt-tl  >H-(E-AINF41,1>  ,E9INT-(t  )) 
DATA  H  T/l  •/ 

— u AT « -  A L 3Hk/  1  »£6  »ci  ■  FOE 7 li  1 » ■«  1  A  3t1tV332 r« 3 ~0 £>■)■§ vrl-S~4»-» 2-? 87E-5) — 
A.766I7E-3/ 


c  INITIALIZATION 


IF ( ( < K  — 1 ) .6E.M)  GO  TO  6 
IF(«.GT.1»  CO  7  0  5 
N1=N0‘-  *2 

— SCAL-C3UM-2/-D&L - 

GO  1  I- 1»  NO 

— uQ -J—J — if  NO - 

PCI, J)-3. 

— COO-H  . .  . . —  — 

A(I,  J)=3. 

--00-1 1-I-l-rNG - 

W(I)  =  C  . 

-C(I)=C  . - 

F  (I)-i,  . 

-P-fl-rl-)  *>»><}  v 0 - 

FU,1)  =  3.C0  2 

-fM2,2F=-J.Ov~2 - 

X3  ( 1 )  •  XTG  (1  ; 

-  -X  4  <  1  >  =  ■)  « - 

tCH(l)=:X3(l  ) 

~EGOM-l->--»(-i - 

it  J-G  • 

•  -«AG=-G-L*  ALFPA<2-J - 

iF(  A=C- .  IT  .-2lv  )  Sl-EXP  ( ARG) 

-  NCIM=N0 - 

CO^O.  (1)  A  LFE  A  <  3) 

-€*<♦=0644  >  ALFhA<»» - 

HID-Umi 

-  IF  ( I  SE  T  .ECU44-W  <-2>-=-Z4 - 

Nl=«+1 

— <2=<K-t - 


X3C<1)  =X3  tKK4*€0b*3frfc - 

IF { ALF  H4( 2)  .EO.j.)  GO  TO  4 

X*»tKlM »4*-X-*f4K K »  4£00*4-tr-S 1  >7  A-LPHAT-24 - 

GO  TO  3 

XM K14-*X4  <  KX  Hfc-OO^O&E - 

CONTINUE 

tCH  ( <K  > * X-34  K  K) — X-*r<4^K» - 

IF(  «.  ij£.  2)  ELGH  <K<)  =  ( EOH  <  KK ) -?  CH (K2) ) 7  CEL 

X30il)«X3<K-l-) - — - 

1FUKK-1)  .  L  £  .K)  GC  TO  150 

•A  <lt-l  - 


* 


A  a,  AlPH4(4) 

A(2,NM»*A(  1,1)  *  AipHA  ( IK  A2 -  - - 

A(2,NLMA(1  2)- ALFt(A(l)*A2 

UC  ?  I»i#hJ2 -  - 

Jl=If2 

- A-C-ii-,-1  >-»RA — - — — - 

2  C0.4Ttt.J-:  -  - - 

CALL  LS IR7  (  KC,b»  CEL,Ef)  ,Eli  INT  »  5  ) 

• --  -  Cfi3=<>-»  ALPhA(U'-A4 - 

bCALl  =  ^:AL-  (  ALP  HA  {  l)-1  A2)  *'  2 

— - 0-C5-(-3t-2-4-*-4  LPRMbU  «-&0hU.-4 - - - 

c 

C  CONPuTE-HEAN -TRA6Cl-N-G-ERROR-(-I-r€-.  -LAG  ANCLE-) - 

C 

. . U=ALPh  A  <  3)-»-H4l  )-f-ALPNA4-4M4H-2)  -f-ALP«A45) 1  4X3(  <K-)-X-4K*H - 

00  lit  1=1,  fC 


0  tI)=C(I)fEAU, 

UC  —  CCNTINUE - 

UC  CCNTINUE 

- f(U=<  l.-CR4-)-XJ-(KtO-fCR4t-X4LKK-) - 

F  (2 J  =  X  3  (KiO ♦  CR3’  (X ♦(«<-«) -X3 (KK-Hl) 

CO  13t  1=1,  NC 
CO  14C  J=l,  i 

- -  -0  ( I  )  =  0  (  I)  f t A INT-  (I,4>*F-(J  J - 

140  CCNTINUE 

- H  ( I  MG  (44 - 

0(I)=G  . 

130 - CONTINUE - 

C 

C-  CGMPU-TS— COV AR-i  ANCE— NATRIX - 

C 

- CQC  (1,-1)  =  (  ALPHA  U )  -♦  ALP  H  A  ( -74— A  &  S  (  E-0H-(-KK)->-*A  La  H  A  (  8  M-ABS  (4-CQh  («.)  U- — 

1  "SCAL 

1) )-SC ALl 

- C  ALL  - M  ULl  (  €A  IHT C0C»NG-»  N1  ,P4-,  IQ) - 

CALL  MJtT (E(,F,  NC,N1,P2,1C1 

- -  00-221— 1=1, NC- - 

CO  22C  J=i,  t>C 

- P-<£,-J  ( IrJ)f  P  2-(-t-rJ-) - 

22  C  CCNTINUE 

150  CONTINUE - 

£LERR=  A ( 1 ) 

. £L30=S0RT(P41rl->-> - ... 

£LTR=M2) 

- kEWWN . . . . 

END 

-  -  SUBROUTINE- G6SA46UZcRR,4-ZS,)t£LG> - - - 

C  CNHON/j/CO Ui ,UEL,<K,N0,Y1H 2),X30 (2) , £L,E-03, A  ZOO,  M, RA, 

1  /?,  NC  1  ,N02  i  N*  A~,UcL  , U»  ELT P,  A  Z-~R,  ISE-T-rZ  1»Z  2 -  - 

0I1CNSI JN  H(  A)  ,P  ( *» »  4)  ,£>1(4i“)  ,  F2(*4,*,)  , ALPmA(  8)  ,A  (4,4) 

- 1  ,-b(  loi  ,E5  ( 1  C-HrE-f-lN?  ( l  H  >■,  f  A  (-.  ,-4 )-,E-A I-N  *-(■*»  ,t )-,^T •* )  r-CCG-(-4,4->-rt-  ( 4) — - 

2,  X3(li  3  I)  ,X*  (13  .i) ,  EON  (11  00) ,  ECCIHIGGU) 

ECUIVtLENCE  (A(l,l>,r,(l)>,  (CA(i,l),Ea(l)>,(EAINT  (1,1I,EP,INT(1») - 

OATA  HT/l./ 

JAVA-  ALPHA  /  *.'57-,  15,,  15  447,- ,0-1 50  6,,  93041,  .1 563  4E*4  ,  ,6419  2  E  *4  ,,005289- 
A/ 


•iV'V*- w*  "sr&~.  ®awni»™5iW”' 


IFKK.GT.l)  GC  10  5 
-fU=H^-»2 - 


O  -INI-T  I AL-J-Z-A-HQ-N - 

C 

- GG — 1 — !■»  1 1  -N  0  ■  ■  ■ 

00  1  J=1,NQ 

■  -  P<I|JH»-T. - 

CCC(It  J)  =  3. 

1 - A(I,J)-=*. - 

JO  11  1=1, NO 
- HH-Hh - 


J(I)-G  . 

- f-l  H  *  C--« - 

11  P  (I,  I)  =  ]*30  C G 

- P  UT  IT*)  «A-i  --£■  -OOSX-EL-G-h-*-2 - 

P»2,2)  =  J.D0  „2  C0S(ELG)<-*2 

- X-C4-l-F*->dG~«M - 

X4(l)s  0 • 

- £0H<  D-aXJ+t) - 

EODH ( 1  )-0  • 

- S13&-, - 

c 

•  C-COHPU  IS-A+rd-irl-Q  A  &  ■  67  A7gS-~X3-A  NO-XA - 

C 

5  - MU)-V-U  U4 - 

6  IF(I3£T.EQ.1>  W ( 2)=Z2*C0S<ELG) 

- C  CNT INUE- - 

CC=C0S ( SLG) 

AF.G=-L'Eu*  ALPrAt  U-C3 
lFtARG.3T.-2Ci .)  S1*EXP< ARG) 

K1=<K+1 - 

K-=«-l 

X3(Kl>-“X3<i<  FHA-ZC-u*9efc - 

1F{AlFHA(1) .EO.C  .)  GO  TO  A 

- X  -  ( K-14-=  3 1— <  XM-X*M-A-ZOO  '-9  E-b) - 

GC  TO  3 

A  A  *•<  Kl)  =XMKK)*AZGGJOEL - 

3  CONTINUE 

C  COMPUTE  ANJ  STORE  ESTHATE0  TARGET  VELOCITY  AND  ACCELERATION 

C - 

cCH(Kl)=X3(K.  l)-XA(Kl) 

I F  (K2 . GE .  1 J  EGG«i  <KK  >=  ( EDM (Kx  )  «€)H  (K2 )  )  /  CEL - 

X  30  t  2 >  =  <3  ( Kl ) 

IF(  (K»  -l)  ,Li.m--W-  TO  -150 - 

LO  2  1=1, N02 

~  J  1=  I-+*2— ■  ■■ — ■■  - - —  -*-•  ■"--■i  ■  .  ■  ■—  —  .1-—  —  ■  . . . . 

A<Jl,I»=ftA 

- - - - 

1  tiEr-J*  (ELG-OE«.GF/OEt - - 

Cfc=COS.(ELGJ 

- i-B*-rhEGD^TAMEL-<r> - — _ 

oCAL  =  ( 0(2)  CL)  *• 2/OEL 

CCR'CC  (  2)  -  AlFi-h-J  2)  aC8 . . . 


M(ltl)  =-C.'R»lt. 

•  -  4(l,?l»-C-j<  c)-ALPHA<3>*€9 - 

A  »2,N01)  =  -Cifi- A2 

- n-i  *U~h*AH  i - 

m (2* 2) =TB 

C  COMPUTE  TRANSITION  MATRIX  EA  ANJ  ITS  INTEGRAL  El I  NT 

C - 

CALL  0SCRT<NC,^,CEL,EB,£3INT,5> 

- Gfc** GC  (  i)  •  ALfhA  ( «t>-* gg - 

CGC(2,  2)=ALPhA(5)*SCAL*A2*»2 

c  J - 

C  COMPUTE  MEAN  TRACKING  ERROR 

c  - 

U-ALPh  A  (2 ) ;  W(l)  ♦ALPhA(3)*M2H-ALPHA{L)‘  (X3< <K) -X L(<K) > 


uO  12C  J=i,  EC 

-  -  LMI  >  =0  H H-E-AHi  ,-J  T'-W-W) - 

12 G  CONiINUE 

11  G  — CONTINUE- - 

F  (l)-tCT-CRL)  y.T  { KK)  ♦  CRL*  X4( KK ) 

- ,  _  x-3-fi^K-f— €r  ♦  frftv‘H*M-KK---Hr— )r3-(**-M  H-*-  A2 - 

00  13f  1=1,  NC 

--  - 00 — 1  *  L — J -  1 - 

U(I)-C(I)+t  aim  <  I,  j»-»f  (  ji 

1L0 - CONTINUE - 

W  (I)- j(I) 

- - 

13C  CONTINUE 

c  - - 

C  COMPUTE  COVARIANCE  MATRIX 

C - 

CQCU,  l)=(ALFhA'5)«-ALPMA(6)  Ac  SCEDh(KK)  )  «-AL9HA  (7  )  C  AMS  (E  ODH 

- fr- *SO*fc - 

CQC<2, 2)= (ALFHA(5)»AlPHA(6)  AGS  (EDH (KK- M > 

.  1---«-AlFHA<7>— Ar,V<£G0H(KK-M)-H*-SGAL»A2«2 - 

CALL  MUlT(EAIM,CCC,N0,N1,P1,1G) 

. C  AL  L  -  M  Uc-T-  (  6  A-,P  rN  OVU-yMriH - 

CO  22C  l=l»  NC 
00  22 C  J*  1, I^C 
P<It  J)  =P1  (I  t  JJ  *P2(I,J) 

220  CONTINUE - 

IsO  CONTINUE 

OELG-»tL‘> - - - - - 

«ZE  RR=  2  (1 )  7CE 

— - A  <S0»S«HTHP-Ut-1'»  )/00 - 

AZTR=W (2)/Cc 

RETURN - 

END 

-  iUBROUTINc-R-LL-^-l  tfP-»CTC-lTHrMR-) - 

DIMENSION  E 1 1),F  (1) ,G(16) ,H(  1) 

. -  OG— lr> -I-trfc - 

11=1 

GO  1-j  K=1tU - 

TEMP=c . 

„  U0-5-J=ItClt-L - 

T  EMP=T  £  1°  ♦£  (  Jl  VF  (II) 


CK< ) >  > 
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1-1*  tit-1 - 

KK  =  (  <-  1 )  *■  L  F  I 

M «<K)iTcMP - 

G <KK)  =  TEMP 
IF { . £1.  llfEHNM- 
00  23  I  =1 ,  L 
0  0—^  0 — K~~I~  i  L  • 
TEMP=t  . 


UC  13  J=I,L1,L 
TEMP*  TEMP t-G-M-JJ-HH-H- 
II-IIt-L 


H  (<K)  =  T clip 


00  33  1=1,  L 

•  -»  —  r  ■  » 


1/0  3 J  J=L3*  L 

K 1=  <1-1  PM.-* - - — - - - - 

K  2=  (  J- 1)  *  L  *  1 

HK1I^H(«) - 

END 

SU-3R0UTlN5-t£frPT-(N^I}4rA-r9fctr,£A-T£«»*TrNT. 
DIMENSION  A ( 1) ,£ A(l), EAINT(l) , CC£F<3C> 


NCIHlsNDIMM 

-  -  NN=NOI  M  -*NOI-* - 

N  TM1  =  NT -1 

-  CC£F<NT>*-l-» - - - 

00  10  1=1, NTH 

- H«  NT  -  I - 

1G  COEF  ( I  U  =  0£l.^C0rF(II+l  J/FLOA  T  ( I) 

- - Nr-  MUiT— feg— AT-tr£-AST---3- - 

CALL  t I AG ( NC IM, t AI NT , A ,COEF { 1 ) ,C0Er (2)) 

-  00  6)  L~*3»NT - 

OrtLL  MULT  (A,  cAINT * N0IH,NN»£A  ,  1) 

- iFU-.M.  N  I-3-&C-I-0— 7-3  - - 

fcO  CALL  OIAGINCIM,  tAlM,EA,l  .0,COEF(L») 

70  L/C  f  0  I  It  1 ,  NN,wu  IH1 - 

E  A(  £  I )  =iIA  ( 1 1 )  tl ,  it 

60  -CCNTlNUE - - - — — - 


DIMENSION  A  (1)  ,«i  (1) 

- -.—NCI  Nl=  N  GI-M  4-1 - - — 

NN=NOIN‘NOI 

. . NM1— NG-tN-i - - 

11=1 

IF(C1.£0.1«  )  GO  TO  10 
DC  5  J  =  l,  NN,  Mi  IN 

- - Kt-ltNMl - 

UO  «♦  I  =  J,  K 

- <*--A  <I>-=C-l  -*fe-<T-) - 

A  cm  =  A  (XI*  <C2 

- 5 — M—-tI  t-MOI  Ml-  ■  — - 

nETLRN 

1C  -  OC-7  -J=  1,  NN-t-M/tN - 

<=J+NM1 

■  00-  6  -I-s  Jr  K- - 

6  A (I i =t  (  U 

- M  M-II-H-  A  MM  i-*-GS - 

7  I 1=1 I ♦ NJI Ml 

-  -  PETU3N - 

ENO 

22C6,3u,G,  2 - - - - - - 
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